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ABSTRACT 


This report describes an analytical study of the performance of flat, 
hydrogen-cooled heat exchanger panels, as part of a comprehensive inves- 
tigation of problems associated with design and fabrication of efficient 
regenerati vel y cooled structural panels. Heat exchangers were designed 
and performance was calculated for several panel concepts in a range of 
loads from 10 to 500 Btu/sec-ft 2 (114 to 5680 kW/m 2 ). Factors that influ- 
ence heat exchanger and manifold system design, such as state of the art, 
materials, and fabrication techniques are set forth. Parametric data on 
the performance of a broad array of coolant passage geometries are pre- 
sented. The potential for coolant conservation and reduction of tempera- 
ture differentials by use of various methods of flow manifolding, flow 
folding, multiple stacking of heat exchangers, and insulation are discussed. 
Analytical and experimental studies of rectangular and triangular flat- 
manifold systems show the effects of geometry variation on pressure drop 
and panel flow distribution. 




FOREWORD 


This report was prepared by AiResearch Manufacturing Company, a division 
of The Garrett Corporation, Los Angeles, California, for the Langley Research 
Center of the National Aeronautics and Space Administration. Results of an 
analytical study performed under Task Order No. 2, Heat Transfer and Fluid Flow 
Analysis of Heat Exchanger Surfaces and Manifold Systems, are presented as part 
of a comprehensive analytical and experimental study of regenerat i vel y cooled 
panels accomplished under Contract NAS 1-5002. This program was under the 
cognizance of Dr. M. S. Anderson and Mr. J. L. Shideler of the Aerothermoe last i - 
city Section, and Mr. R. R. Howell and Mr. H. N. Kelly of the 8-Foot High- 
Temperature Structures Tunnel Branch of the Structures Division, Langley 
Research Center. 
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HEAT TRANSFER AND FLUID FLOW ANALYSIS OF 


HYDROGEN-COOLED PANELS AND MANIFOLD SYSTEMS 

By F. M. Walters and 0. A. Buchmann 
AiResearch Manufacturing Company 
a division of The Garrett Corporation 


SUMMARY 


As part of a comprehensive investigation of problems associated with 
design and fabrication of efficient regenerat i vel y cooled structural panels, an 
analysis was made of the heat transfer and the fluid flow performance of flat, 
hydrogen-cooled heat exchangers. In addition, the fluid flow performance of 
selected manifold concepts was investigated analytically and experimentally. 

Within ground rules established as compatible with the state of the art for 
materials and fabrication techniques, a large variety of conceptual designs for 
heat exchanger panels were investigated to determine specific design features 
and ranges of applicability. Performance of a broad array of heat exchanger 
geometries was determined over a range of heat fluxes to 500 Btu sec-ft 2 
(5680 kW/m 2 ), and is presented in terms of coolant flow rate, flow length, cross- 
sectional temperature difference, and pressure drop. Potential advantages of 
various heat exchanger geometries and methods of flow manifolding, flow folding, 
multiple stacking of heat exchanger elements, and insulation were investigated 
qualitatively and quant i tat i ve 1 y, in terms of coolant conservation and tempera- 
ture differential reductions. 

In general, the investigations were made assuming fixed uniform heat fluxes 
over the exposed surface of the panel. However, for coolant conservation 
studies a uniform recovery temperature was assumed. A hydrogen coolant inlet 
temperature of I00°R (55.5°K) and various outlet temperatures from I400°R to 
I900°R (778° to I055°K) were used. The inlet pressures were varied, as 
requi red--subject to a maximum pressure limitation of 1000 ps i (6880 kN/m 2 )--to 
maintain a coolant outlet pressure of 250 psi (1720 kN/m 2 ). 

Although the single-pass heat exchanger was found generally applicable to 
the entire range of conditions studied, two heat exchanger concepts appeared 
attractive for coolant conservation: (l) insulation on the heat exchanger hot 
surface and (2) flow folded in the panel width dimension with counterflow in 
alternate passages. Typical hydrogen flow rate reductions of 5 to 10 percent 
are achieved by use of folded flow at heat fluxes below 100 Btu/sec-ft 2 (1135 
kW/m 2 ) and recovery temperatures below 7000°R (3890°K). Hydrogen flow reduc- 
tions of 20 percent or more can be achieved by use of insulation throughout the 
range of heat flux for this study. 



In addition to the parametric analysis of heat exchanger panels, detailed 
performance was calculated for selected panel concepts. The concepts evaluated 
were (l) a single-sandwich configuration that provides for both structural load 
carrying capability and coolant containment and flow routing, (2) a composite 
configuration in which the coolant-pressure-containing surface heat exchanger 
is metal lurgical ly bonded to the structural-load-carrying sandwich panel, and 
(3) a coo 1 ed- sh i ng 1 e configuration in which the surface heat exchanger is mechan- 
ically attached to 1 ow-tempe rature load-carrying structure. 

The manifold analyses and testing performed as part of this program 
demonstrated that serious flow distribution problems can occur in regenerat i vel y 
cooled panel applications. Geometry variations affecting manifold pressure drop 
and heat exchanger panel flow distribution were investigated experimentally, and 
means of improving the flow distribution are indicated. 


INTRODUCTION 


Vehicles traveling at hypersonic speeds are exposed to severe thermal 
environments. Maintaining the temperature of structural elements of these 
vehicles within the limit of current material technology becomes a serious prob- 
lem. For spacecraft and research aircraft where the heating period is short and 
refurbishment is not a problem, acceptable thermal limits have been maintained 
by designing the vehicle as a heat sink or by ablatively cooling the surface. 

For hypersonic cruise aircraft where heating periods are much longer and refur- 
bishment is a problem, these means of thermal protection are no longer 
sufficient. 

Major portions of a hypersonic aircraft can and will be radiatively cooled; 
however, in many areas the radiation equilibrium temperatures will exceed the 
material limitations and some active means of thermal protection must be pro- 
vided. One of the more attractive means of active cooling for th i s app 1 i cat i on 
is regenerative cooling (use of the fuel as a coolant). This technique appears 
especially attractive for the hypersonic cruise vehicle because cryogenic hydro- 
gen, which has been proposed as a fuel, is an excellent coolant. 

Regenerative cooling has been used successfully for hydrogen-fueled rocket 
engines. However, in contrast to rocket engine applications, which are charac- 
terized by small areas, high heat fluxes, and short operating times, air breath- 
ing hypersonic cruise aircraft will have large areas of low to moderate heat 
flux and will be expected to operate for much longer periods of time. Conserva- 
tion of coolant and minimization of weight become paramount for these vehicles. 

Therefore, a study has been initiated to investigate problems associated 
with design and fabrication of efficient regenerat i vel y cooled structural panels 
The overall study is concerned with practical engineering problems of material 
applications and fabrication, structural design and analysis, heat transfer and 
fluid flow analysis, and the integration and interaction of various aspects of 
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the general design problem. In an initial investigation, reported in ref. I, 
analytical studies were conducted of a wide array of conceptual designs for 
flat, hydrogen-cooled panels operating at heat fluxes from 10 to 500 Btu/sec- 
ft 2 (114 to 5680 kW/m 2 ) and normal pressures from 6.95 to 250 psi (48 to 1720 
kN/m ). Procedures for integrating the heat transfer and structural design and 
minimizing configuration weight were developed, and ranges of applicability of 
the various conceptual designs were determined. 

In a companion investigation, reported herein, the heat transfer and fluid 
flow performance of flat, hydrogen-cooled heat exchanger panels was studied 
analytically; the performance of associated manifolding systems was studied 
both analytically and experimentally. Various influences on heat exchanger 
performance were examined and potential advantages of various methods of flow 
manifolding, flow folding, multiple stacking of heat exchanger elements, and 
insulations were investigated qualitatively and quantitatively, in terms of 
coolant requirements and heat exchanger temperature differentials. In addition, 
parameteric data were generated on the performance of circular, rectangular, 
interrupted rectangular, and pin-fin coolant passages. 

Although the present investigation was closely allied with the configura- 
tions of ref. I, the types of problems encountered and investigated are 
representative of those that will arise in the design of hydrogen-cooled panels 
for application to hypersonic vehicle or engine surfaces. Furthermore, the 
data presented encompass the ranges of app 1 i cab i 1 i ty for both exterior and 
interior surfaces and provide the basic quantities of interest for specific 
design or tradeoff studies. 


SYMBOLS AND PARAMETERS 

A = area, ft 2 (m 2 ) 

A = fluid designation in four-fluid heat exchanger program 

b = fin or web spacing, or width, in. (cm) 

C = ratio of heat exchanger thermal conductances at 2000°R (|||| 0 K) and at 

T co , dimensionless; fluid designation in four-fluid heat exchanger pro- 
gram 

C p = specific heat at constant pressure, Btu/lb-°R (j/g-°K) 

D = diameter, in. (cm) 

f = fin effective thickness factor, or friction factor, dimensionless 

G = fluid designation in four-fluid heat exchanger program 

9 = 9 C = 32 * 2 lb (mass) ft/lb (force) sec 2 (980.7 cm/sec 2 ) 
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H 

H 

h 

j 

K 

k 

L 

o 

1 

N 

P 

Pr 

Q. 

q 

R 

Re 


enthalpy, Btu/lb (J/g) 

fluid designation in four-fluid heat exchanger program 

height, in. (cm) or heat transfer coefficient, Btu/sec-°R-f t" 
( kW/° F-m 2 ) 


Colburn modulus, St(Pr) z/3 , dimensionless 

pressure drop coefficient for turning, expansion, or contraction 

thermal conductivity, Btu/hr-°R-ft (W/m-°K) 

offset length (uninterrupted flow length), in. (cm) 

panel length, or hydrogen flow length, in. (cm) 

effective fin length, in. (cm) 

number of fins/unit of heat exchanger width 

pressure, psi (kN/m 2 ) 

Prandtl number, |j,Cp/k, dimensionless 
velocity head, W 2 /A 2 ^.2g^_p, lb/in 2 (kN/m 2 ) 
heat transfer rate, Btu/sec ( kw) 

resistance to flow, dimensionless; gas constant in equation for 
dens i ty 


Reynolds number, 


4 r, W 
h 

mV 


dimensionless 


r 


St 


= radius, in. (cm) 

= hydraulic radius, in. (cm) 


= Stanton number, 



dimensionless 


T = temperature, °R (°K) 

TC, TH = fluid temperatures in four-fluid heat exchanger program, °R (°K) 

T I T2 

, , \ = metal temperatures at nodes I, 2, ..., in four-fluid heat exchanger 

program, °R (°K) 

t = thickness, in. (cm) 


4 


W = flow rate, lb/sec (kg/sec) 

WC, WH - fluid flow rates in four-fluid heat exchanger program, Ib/sec (kg/sec) 
w = panel width, or hydrogen flow width, in. (cm) 

X = variable width, in. (cm), in dimensionless core flow width, X/W 

a = thermal expansion coefficient, in./in.-°F (m/m-°K) 

y = material density, lb/in 3 (kg/m 3 ) 

A = change in, increment of 

e = strain, in. /in. (m/m) 

Tj = fin effectiveness 

1\ = overall fin effectiveness 

p, = viscosity, Ib/sec-ft (kg/sec-m) 

p = density, lb/ft 3 (kg/m 3 ) 

p = density at average pressure and temperature 

a 

a = stress, psi (kN/m 2 ), or ratio of density at average temperature and 

pressure, lb/ft 3 to standard density of I lb/ft 3 

Subscri pts : 

a = average 

C = coolant 
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HF 

HM 

HT 

I 

i = 

LP 

M 

m = 

max = 

min = 

0 

o = 

P 

R 

req = 

SP 

std = 

T 

W * 

WH 

i > 3^ 

4 , 5 = 


horizontal fin 

hydrogen temperature between folded and single-pass panels 
hydrogen temperature at |80°-turn in folded-flow routings 
inlet 

insulation 
longest path 
manifold 
momentum 
maximum 
minimum 
out 1 et 
overa 1 1 
port 

recovery or adiabatic wall 

requ i red 

shortest path 

standard 

total 

hydrogen-cooled surface of hot wall 
aerodynamic-heated surface of hot wall 

flow distribution test unit test station numbers 
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Fin geometry is designated with a six-part nomenclature: 
20 (7 . 9 )R- . 1 00 (. 2 54 )-. 1 25 ( . 3 I 8 )- . 004 (. 010 ) ( 10-61 ) 


1 Figure number from 

reference 2 used for 
performance calculation 

l t f 

in, in. (cm), fin thickness, 
tube wall thickness or fin 
d i ameter 

Lqj (cm), uninterrupted flow length, 

full wave length for wavy fin, L Q /4 r^ = 

100 for tubes and plain fins, staggered 
row center spacing for pin fins 
^fin* ( cm ); height or tube outside diameter 
Designation of cross section and/or type of flow length interruption 

R = rectangu 1 ar fins 
T = t r i angu 1 ar fins 
W = wavy f i ns 
TB = tubes 
P = p i n fins 


L Fi n spacing N = l/b f . , fins/in. 
of tubes per in. (cm) 



U- b fin-J 



( f i ns/cm)( for tubes, N = twice the number 



HEAT EXCHANGER GEOMETRY NOMENCLATURE 
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STATEMENT OF PROBLEM 


Objecti ves 

The objectives of this study were to examine the heat transfer and fluid 
flow characteristics of heat exchangers and manifolds for flat, hydrogen-cooled 
panels and to provide useful data for the design of panels of this type. The 
desired capabilities of an efficient heat exchanger panel and manifold system 
were considered to be: (l) to limit maximum temperature for strength or cor- 

ros ion resistance, ( 2 ) to limit steady-state AT th roughou t the structure f o r 
acceptable thermal stresses, (3) to limit coolant flow rate and, thereby, 
increase efficiency and reduce vehicle weight, (4) to make efficient use of pres 
sure drop for obtaining maximum heat transfer coefficient and limiting required 
pumping power, and (5) to establish flow distribution to fit the heat flux dis- 
tribution across the panel width. 


Environmental Conditions and Design Constraints 

During this study, the environmental conditions applied were considered 
representative of conditions that may be encountered on external and internal 
(i.e., inlet, duct, and engine wall) surfaces of hypersonic cruise vehicles. 

The design constraints imposed were those considered representative of good ^ 
engineering practice, and they were based on present day materials and fabrica- 
tion technology. Limiting conditions used in the study presented herein are 

• Panel length--up to 5 ft (1.525 m) 

• Panel conf i gu ra t i on-- f 1 at 

• Coolant--hydrogen 

• Cooling method — forced convection 

• Net heat flux — 0 to 500 Btu/sec-ft 2 (0 to 5680 kW/m ) 

• Coolant pressure--250 to 1400 psi (1720 to 9650 kN/m ) 

• Coolant temperature-- 100° to I900°R (55.5 to 1055 K) 

Heatinq condition.-- For the purpose of this study, it was assumed that 
heating occurred from one side only. When two-sided panel heating occurs, heat 
fluxes will generally differ significantly on the two sides and special coolant 
controls will be required to permit two-sided heating of a single panel with 
acceptable thermal stresses. Also, installation requirements may lead to the 
use of separate panels for each of two surfaces with different heat fluxes. 

Two different heating conditions were considered. For the majority of 
cases, a uniform heat flux was assumed over the surface of the panel. ^ Due to 
temperature variation along the length of a cooled panel, this condition is 
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approached for a uniform external environment only as recovery temperature 
becomes very large relative to wall temperature. Therefore, the term infinite 
recovery temperature is used herein to refer to the uniform heat flux condition. 
For the remaining cases, the panel was assumed to be exposed to a hot gas with 
a uniform finite recovery temperature so that the effects of panel surface tem- 
perature could be considered. The nominal level of heating for these cases was 
defined by the heat flux at the cold end of the panel, which was assumed to be 
operating at a hot wall temperature of 500°R (278°K). The actual heat flux 
along the panel varied directly with the difference between the local hot wall 
temperature and the hot gas recovery temperature; therefore, the average heat 
flux to the panel was less than the nominal value. Uniform finite recovery 
temperatures of 3000°, 5000°, and 7000°R (1670°, 2780°, and 3890°K) were used 
for these studies. 

Coolant Pressure.-- For this study, the minimum coolant pressure that 
occurs at the outlet manifold was chosen to allow for a pressure differential 
that can be used to inject the hydrogen into an engine combustor section. The 
use of supercritical outlet pressures allows the assumption of forced-convec- 
tion, single-phase, heat transfer coefficients throughout the study. The 
coolant inlet pressure selected was that pressure required to produce the 
necessary coolant flow through the heat exchanger. The upper limit on inlet 
pressure was usually taken to be 1000 psi (6895 kN/m 2 ), although inlet pressures 
up to 1400 psi (9650 kN/m 2 ) were considered. 

Temperature. --A hydrogen inlet temperature of I00°R (55.5°K) was used during 
this study. It was assumed that hydrogen would be stored at 40°R (22.2°K) or 
less and would undergo a temperature increase due to heat leak and/or compressor 
energy input of about 60°R (33.3°K). A hydrogen outlet temperature of I600°R 
( 889 K) was usually assumed, although outlet temperatures from 1400° to I900°R 
(778° to I055°K) were considered. 

Practical structural operating temperatures have been limited to approxi- 
mately 2000 R (Mil 0 K) by current state of the art in materials and fabrication 
of heat exchangers and structures. This temperature was used throughout the 
program as the maximum permissible temperature at any point on regene rat i vel y 
cooled structures. 

Materials,-- Use of superalloys such as Hastelloy X, Inconel 625, and 
Waspaloy was dictated by st rength-at-temperature and 1 ow-cycl e- fa t i gue charac- 
teristics. The important material property for heat transfer analysis is ther- 
mal conductivity. Thermal conductivities for these three materials are shown 
i n f i g. I . 

Heat exchanger geometry .-- Table I specifies the limitinq dimension for 
each of the heat exchanger geometry variables and the primary limiting condi- 
tions. The table of Symbols and Parameters defines these geometry variables 
and descriptive terminology. 

Hydrogen properties .-- The transport and thermodynamic properties 
referenced throughout this study are those of parahydrogen gas at 500 psi (3450 
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), as shown in fig. 2. The use of parahydrogen properties is based on the 

...... _ I ! ^Iar>f -h i n f C3t3lytlC PGSC" 


kN/m 
ass 

tion between 


umption of liquid hydrogen storage and insufficient time or c ^ a Y 1 
n between storage and panels for conversion to equilibrium hydrogen. 


RESULTS AND DISCUSSION 


Heat Exchanger Performance 

The important variables and parameters related by the various performance 
analyses are 

• Heat flux, q/A, and recovery temperature, Tr 

• Hydrogen flow length, l 

• Hydrogen flow rate per width of panel, W/w 

• Hydrogen inlet pressure, Pqj 

• Hydrogen pressure drop, &Pq 

• Maximum heat exchanger temperatures, TyR and Ty 


• Maximum structure temperature, Tp 

As detailed in Appendix A, operating characteristics of typical panel heat 
exchangers permit simplifying combinations of the seven performa ^es^nd 5 
listed above. However, the important relationships among these variables and 
parameters are well summarized for all flow routings and geometries by the 
following statement: panel heat exchanger performance is good when smal tem 
perature differences result with low pressure drop and low coo n 

R prove rv temoerature effects. -From the basic equation for convective heat 
transfer from a hot gas to an exposed surface (eq. I of Appendix A) it can be 
seen thlt the heat flux to a panel surface is directly proportional to the dif- 
ference between recovery temperature of the hot gas, T R , an t e tempera ure 

the hot surface of the panel, T^. If the hot surface temperature can be 

increased to approach the recovery temperature more closely, the heat: flux to 
the panel will be reduced and consequently (see eq. 2 of Appendix A) the coolant 

flow rate can be reduced. 

Usinq eqs. I and 2 of Appendix A, a cooling efficiency factor has been 
established to provide an insight into the effects of hot gas recovery tempera- 
ture on the thermal performance of cooled panels and to provide a enter on to 
assist if Liluation'of various cooling configurations. The cool , ng efficiency 
factor is defined as the ratio of coolant flow rate required for a pan 
uniform hot wall temperature of 2000°R (IHI°K) to the flow r *‘' ™ d | 0r 

panel with some other average hot wall temperature, T MH . (A temperature of 
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2000°R (llll°K) was selected as a reference because it represents a practical 
maximum temperature for structural applications of superalloy materials.) Since 
coolant flow rate and average best flux are linearly proportional, the cooling 
efficiency factor is also the ratio of average heat fluxes into the two panels. 


The cooling efficiency factor has been evaluated for selected average hot 
wall temperatures and the results are presented in fig. 3 as functions of the 
hot gas recovery temperature. Hot wall temperatures of less than 2000°R ( I II I °K) 
are representative of typical uninsulated panels in which the hot wall tempera- 
ture varies from some low value to the maximum allowable temperature of 2000°R 
(II I I K ) 35 the coolant temperature increases along the panel. An average tem- 
perature of 1000 R ( 555° K) can be considered as a limiting case for single-pass, 
straight-through cooled panels where the hot wall surface temperature varies 
linearly from 0°R (0°K) to 2000°R (llli 0 K). Average temperatures between I000°R 
(555 K) and 2000 R ( I I I I 0 K) can be obtained with sophisticated flow routings or 
coolant passage contouring, which permits the surface temperature to approach 
the permissible maximum nearer the inlet end of the panel, or through the use of 
insulation on the hot face of the panel. Average temperatures greater than 
2000 R (Mil K) cannot be obtained with flow routing or passage contouring alone; 
these temperatures are representative of cases where nons t ructura 1 insulating 
material has been added to the hot face of the panel to increase the surface tem- 
perature. A temperature of 2500°R (I389°K) appears to be a practical maximum 
temperature limit for nonst ructural metallic type insulation for design life in 
the tens of hours. 


It is apparent from fig. 3 that the cooling efficiency is a relatively 
strong function of wall temperature at the lower recovery temperatures, but 
becomes insensitive to wall temperature at the higher recovery temperatures and 
converges to a cooling efficiency of 1.0 at an infinite recovery temperature. 

As a consequence, an infinite recovery temperature has been assumed for a major- 
i ty of the heat transfer calculations because it permits the heating conditions 
to be defined by heat flux alone and avoids the iterative, incremental analysis 
required to account for variations in the hot surface temperature along the 
panel length. At recovery temperatures of 7000°R (3889°K) or greater, the cool- 
ant flow rate obtained by assuming the entire surface of a panel is at the maxi- 
mum allowable temperature of 2000°R ( I I I I 0 K) differs by less than 15 percent 
from that obtained when surface temperature variations are considered. For 
lower recovery temperatures, the differences resulting from such an assumption 
become progressively larger and more detailed incremental calculations involving 
the hot wall and recovery temperatures are required for accurate prediction of 
coolant requirements. 

The recovery temperature of the hot gas must always be specified to evalu- 
ate the coolant conservation potential of various flow routings and insulations. 
It is apparent from fig. 3 that the potential for improvements in cooling effi- 
ciency due to sophisticated flow routing or insulation decreases rapidly as the 
recovery temperature Increases. For example, at a recovery temperature of 5000°R 
( 2778°K) , an increase in the average surface temperature from I000°R (555°K) to 
3000°R (|667°K) will increase cooling efficiency by approximately 75 percent; 
whereas, at a recovery temperature of 7000°R (3889°K) the same change in surface 
temperature nets a 45-percent increase in cooling efficiency. It should be 



noted, however, that the quantity of coolant conserved is the same at both 
recovery temperatures. Thus, although the percentage of coolant that can be con- 
served by a given increase in average hot wall temperature varies with recovery 
temperatures, the actual quantity of coolant that can be conserved is direct y 
related to the increase in hot wall temperature that the flow routing or insu- 
lation provides and is independent of recovery temperature. Furthermore, or a 
qiven hot wall temperature, the heat flux to the panel increases as the recovery 
temperature increases; consequently the thickness of insulation requi red to 
attain that hot wall temperature decreases. Therefore, although the coolant 
efficiency decreases, the efficiency of the insulation in terms of coolant reduc 
tion per unit weight of insulation actually increases with recovery temperature. 


Single -pass heat exchanger panels .-- As discussed in detail in ref. I, 
panel stress is caused by normal pressure, dynamic pressure, coolant internal 
pressure, and various temperature gradients, both transient and steady-state, 
which may exist in all three panel dimensions. The temperature gradient through 
the heat exchanger thickness causes important thermal stresses and is defined as 
the sum of the temperature differences through the hot face sheet or the tu e 
wall, T uh - T j and the fins, T - T , as indicated in fig. 4. 


For a 0.010-in. (0.025-cm) thick superalloy heat exchanger face sheet, 

T - T ranges from I60°R (89°K) at the cold end to 100 R (55.5 K) at the hot 

end with a heat flux of 500 Btu/sec-ft 2 (5680 kW/m 2 ). The face sheet AT varia- 
tions are due to superalloy thermal conductivity variations with temperature, 
and are obtained as specified in Appendix A. 

In ordinary heat exchangers, where the heat flux is lower tt ?^ t 

considered here, the fin AT is much smaller than the wal 1-to-bul k-f 1 uid AT. In 
fact the fin tip most remote from the heated surface is at approximately the 
same temperature as the directly heated prime surface. However, for heat flux 
rates above 10 Btu/sec-ft 2 (114 kW/m 2 ), the fin tip attached to the back sheet 
is essentially at the coolant local bulk temperature and T y - T p is equal to or 

greater than 95 percent of - T c . At heat flux near the high end of the 

design range for this program, approximately one-half of the entire fin is at 
the local coolant temperature. Therefore, since it was assumed that there was 
no heat loss out of the backside of the structure supporting the heat exchanger, 
the temperature of this structure was the same as the local coolant temperature 
and the maximum hydrogen temperature was also the maximum temperature of the sup- 
porting structure. The primary thermal stresses occur because the hot face sheet 
of the heat exchanger, which i s at a higher temperature, is constrained to remain 
flat and hold the planform dimensions of the supporting structure. 


The basic function of heat exchanger geometry is control of thermal 
conductance and pressure drop. Thermal conductance is that property of the 
cool ant-and-geomet ry combination that controls temperature differential between 
the hot-gas heated surface and the structure, i.e., T^ - T p . Pressure rop an 

outlet pressure together control pressure-containment strength requirements 
and determine pumping power required for panel systems. The heat exchanger 
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geometry must provide a lightweight and fabricable configuration, contain the 
hydrogen pressure, maintain acceptable heated-surface temperatures, and provide 
cross section temperature differences within the thermal stress limitations of 
available materials. 


The heat exchanger geometries considered for use in hydrogen cooled panels 
are shown in fig. 5. Geometries with discrete passages such as round and rec- 
tangular tubes may be distinguished from those geometries formed from corrugated 
tension members brazed between two flat sheets. Smooth aerodynamic surfaces 
can be obtained with the heat exchanger geometry in which corrugated metal fins 
are brazed between flat sheets. The aerodynamic and heat transfer effects of 
the rough surface presented by tubes are adverse factors in some applications. 

The important geometry characteristics that allow variations in thermal 
conductance and pressure drop, and which apply equally to all geometries, are; 
(I) cross section shape, (2) uninterrupted flow length (offset length), (3) 
types of flow length interruption, ( 4 ) fin height, (5) fin spacing, (6) fin 
thickness, and ( 7 ) thermal conductivity. "Tube" or "tube wall" can be substi- 
tuted for the word "fin" in the previous sentence. 


Because heat exchanger geometry efficiency is measured in terms of 
hydrogen flow rate, pressure drop, and material weight, high thermal conduc- 
tance at a given pressure drop is usually sought. High thermal conductance 
reduces fin AT and permits the hydrogen temperature to be maximized within the 

constraints of the maximum heat exchanger face-sheet temperature, T , and the 

WH 


maximum structural 


tempe ratu re, 


T 


F* 


Fin AT is reduced and pressure drop is increased ( I ) by reducing 
un i n te r rupted flow length, fin height, and fin spacing, or (2) by Increasing 
fin thickness. Frequent interruption of the hydrogen boundary layer along the 
passage walls (e.g., by use of offset fins) can be used to achieve high thermal 
conductance. Increasing fin thermal conductivity will decrease fin AT without 
increasing pressure drop. Although discrete passage geometries inherently can 
withstand internal pressure better than corrugated tension members brazed 
between two flat sheets, the typical hydrogen inlet pressures required for 
regenerative panel applications are within the pressure capabilities of the 
latter conf i gurat i on. 

For the tubular heat exchanger geometries without braze fillets, the 
corrugated surface area is 57 percent greater than the planform area, w£ . In 
the limit, an increase in heat flux and hydrogen flow rate of 57 percent over 
flat surfaces can be postulated. In ha 1 f- tube-d iameter braze fillets, as 
assumed, and as defined in the Symbols and Parameters table, exposed area is 
18 percent greater than w &. In this case, heat flux and hydrogen flow rate 
increases of up to 18 percent greater than those for flat surfaces can be postu- 
lated. Hydrogen thermal conductance was estimated for both of the above surface 
condition assumptions; however, no increase in heat flux or hydrogen flow rate 
was assumed because of increased exposed area. Hydrogen pressure drop is not 
affected since plain, round, inside surfaces were assumed for all tubes. Tubes 
have higher thermal conductance without braze fillets than with braze fillets. 
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Without braze fillets, half of the tube inside circumference was assumed to be 
fin area transferring heat to the hydrogen. With the more realistic assumption 
of half-tube-diameter braze fillets, three-fourths of the tube inside circumfer- 
ence was assumed to be fin area transferring heat to the hydrogen. The fin AT 
values, T - T , were assumed to occur in one-fourth of the tube inside circum- 
ference without braze fillets and in three-eighths of the tube inside circumfer- 
ence with half-tube-diameter braze fillets. By compa r i son, f i n length for 
rectangular fins is taken as the fin height plus half the fin spacing. 

The above considerations led to emphasis on the rectangular-offset-fin 
configuration throughout the analyses. Work was done with the other geometries, 
however, to point out their strengths and weaknesses, but they were not studied 
extens i ve 1 y . 

The performance of heat exchangers with single-pass flow routing was 
calculated by the method described in Appendix A for the seven important geo- 
metry variables previously listed. A somewhat smaller range of geometry 
variables was examined more intensively (Appendix B) to provide temperature and 
pressure data for heat exchangers studied in ref. I. Figs. 6 through 15 show 
heat exchanger performance in terms of hydrogen thermal conductance per unit of 
panel projected area and density-adjusted pressure drop perunit of hydrogen 
flow length, both as functions of hydrogen flow rate per unit of panel width. 
Plots of density-adjusted pressure drop and fin AT as functions of heat flux, 
with panel length as a parameter, are presented in figs. 16 through 20. Fig. 

21 shows the effect of temperature-dependent fluid properties on thermal conduc- 
tance and density-adjusted pressure drop. 

Since only friction pressure drop is presented in the performance curves, 
overall panel pressure drop must be calculated by adding man i fold pressure drop 
as well as heat exchanger core entrance, exit, and flow-acceleration pressure 
drops. The flow-acceleration pressure drop in hydrogen-cooled panels is about 
6 percent of the friction pressure drop. The core entrance and exit pressure 
drops are usually so small as to be negligible for the panel geometries con- 
sidered. They will, of course, represent a somewhat larger percentage of the 
friction pressure drop for smooth tubular surfaces than for frequently inter- 
rupted surfaces such as rectangular offset fins. To make the heat exchanger 
pressure drop curves applicable at any hydrogen gas density, all pressure drops 
were multiplied by a density ratio. Viscosity variations also affect frictional 
pressure drop; however, as indicated by the small variations wi th. temperature 
of the density-adjusted pressure drops for offset fins shown in Figure 21, the 
effects of viscosity variations are insignificant for the temperature range and 
flow rates of interest. At flow rates lower than those shown in fig. 21, heat 
exchanger pressure drop is less than 2 percent of a typical panel inlet pressure 
of 500 psia (3450 kN/m 2 ). 

Flow rout i ng . --Cool i ng flow rate can be reduced (cooling efficiency 
improved) by increasing average wall temperature, T^. Use of a more complex 

hydrogen flow routing than the single-pass concept (fig. 22a with one inlet) 
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produces higher wall temperatures and higher cooling efficiencies. The single- 
pass flow routing concept is distinguished by: 

• Coolant flow parallel or counter to external hot-gas flow 

• Hydrogen flow length equal to panel length 

• Inlet and outlet manifolds located at opposite ends of the panel 

• Average wall temperatures between 1000° and I500°R (555° and 
834°K), depending upon hydrogen outlet temperature and heat flux, 
under the boundary conditions specified for this program 

Preceding a detailed discussion of the flow routing concepts depicted in 
fig. 22, some of the results of flow routing studies are presented to summarize 
cooling efficiency. Table 2 shows that flow routings more complex than single- 
pass yield significant increases in cooling efficiency only at low flux and low 
recovery temperature. At and above recovery temperatures of 5000°R (2780°K) and 
heat flux of 250 Btu/sec-ft 2 (2840 kW/m 2 ), cooling efficiency is essentially 
equal for all flow routings, although the folded- in-width configuration (fig. 
22d) is best. The cooling efficiencies at recovery temperatures of 3000° and 
7000 R (1670 and 3890°K) resulted from preliminary analyses of heat exchangers 
W L t c rt ^o qU fl f1 ° w length. The cooling efficiencies at a recovery temperature 
of 5000 R (2780 K) result from more exact analyses of heat exchangers with 
equal flow length. An example illustrating the method of heat exchanger system 
design for folded flow, as used to obtain the cooling efficiencies at recovery 
temperature of 5000°R (2780°K), is presented later in this section. Although 
flow folding increases cooling efficiency over only a limited range of the 
boundary conditions that were investigated in this study, its potential use- 
fulness will depend upon the following considerations: 

• The extent of the vehicle or engine surfaces operating at low heat 
fluxes. In general, large portions of external surfaces and engine 
Inlets will be operating at low heat fluxes. Given this situation 
and a requirement for regenerative cooling, significant reductions 
in coolant flow requirements could be achieved by flow folding. 

• The complexity involved in providing for folded flow, 

• Off-design performance of and coolant control in folded-flow panels. 

Multiple Inlets: Heat exchanger length in a single panel may be divided 
with all parts connected in parallel^ as shown in fig. 23^ or with parts con- 
nected in series,, as shown in fig. 22a. Compared to single-pass flow, connect- 
ing parts in parallel provides a large pressure drop reduction and no Increase 
in cooling efficiency. Series connection of the parts produces significant 
increases in cooling efficiency and some reduction in pressure drop. 

Where the flow length is a fraction of panel length, the number of inlet 
and outlet manifold sets may be equal to the number of flow length subdivisions 
as shown in fig. 23a. The total number of manifolds can be held to a minimum 
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by using common inlets and outlets as shown in fig. 23b. While thermal stresses 
are less severe with common inlets and outlets, the only structura y accep 
configuration generally has an inlet at each end of the panel and a common hot 
outlet at the center of the panel. In addition, the panel 1 eng th-to-wi t 
ratio should be at least 2. Since a panel length of 3 ft (0.915 m) and pressure 
drop of < 700 psi (4840 kN/m 2 ) at P CQ = 300 psia (2070 kN/m ) can be obtained 

with the single-pass configuration at the highest heat flux considered in this 
study, the need for manifolding in the middle of panels is limited to systems 
with lower pressure drop limits. Some advantages in reduced manifold weight 
may be obtained by use of midpanel hot manifolds. 


Increasing average wall temperature by injecting part of the total 
hydrogen flow at the inlet and at subsequent points along the length of the 
panel, as shown in fig. 22a, causes the sawtooth temperature profiles depicted 
in fig. 24. The thermal stresses may become unacceptable, but are less severe 
than for fig. 23a. Separating a mul t i pi e- i n 1 et panel into several separate 
single-pass panels does increase efficiency with acceptable thermal stress, but 
is complex. In the limit, an infinite number of injection points or short 
panels will allow the entire surface to be at the maximum temperature, which for 
this study is 2000°R ( 1 I I I 0 K ) . 


Figure 24a depicts periodic injection and removal of hydrogen to keep mass 
flow rate constant, thereby maximizing coolant temperature while reducing 
pressure drop. Increasing fin height in a series of steps (fig. 24b) reduces 
pressure drop by providing greater flow area. For concepts shown in figs. 4a 
and 24b, a temperature sensor located at the end of each section controls a 
cold coolant valve at the beginning of that section. These concepts were rejec 
ted for purposes of this study because of control and manifold intricacy, which 
result in high weight and fabrication complexity. The heat exchanger concept 
of fig. 24c has multiple coolant inlet from a cold heat-exchanger layer that is 
insulated from the hot heat-exchanger layer. Instead of controllable valves, 
the fixed orifice position must be designed to provide the necessary flow rates 
The practical problems of connecting two heat exchangers that are operating at 
very different temperatures make this concept appear impractical for inlets 
spaced at intervals of a few inches. Difficulties in obtaining a match between 
design point flow apportionment and off-design conditions also make this system 
appear unattractive. Again, excessive weight and fabrication complexity caused 
rejection of this concept for purposes of this study. 


Flow folded in depth: The fol ded- i n-depth concept (fig. 22b) increases 

cooling efficiency when heat flux from the hotter layer to the colder layer 
exceeds the aerodynamic heat flux to the hotter layer on a substantial fraction 
of the panel length. By changing heat exchanger geometry, the heat exchanger 
thermal conductance can be changed relative to any given hot-gas thermal conduc- 
tance. Reducing the thermal conductance of the cooler heat exchanger layer 
reduces the hea? transfer from the hotter to the colder layer and thus reduces 
average wall temperature and cooling efficiency. In the limit, the thermal con- 
ductance of the cooler heat exchanger layer can be zero. Then, no heat is tran 
ferred from the hotter layer to the colder layer, single-pass temperature pr 
result throughout the hotter layer, and the colder layer is entirely at the 
hydrogen inlet temperature. 
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For analyses in this study, the sheet metal layer between the hotter and 
colder heat exchanger layers was assumed to have zero resistance perpendicular 
to the thickness. As the resistance of this layer is increased, or if insula- 
tion is used between the two heat exchanger layers, the heat transferred from 
the hotter layer to the colder layer is reduced. Effectively zero resistance 
can be obtained if the hotter and colder heat exchanger layers are connected by 
a single sheet metal thickness. In this case, however, thermal stresses result- 
ing from hydrogen inlet and outlet temperature differentials become a limitation. 
If the hotter and colder fin layers are not metal 1 urgi cal 1 y bonded to a single 
intermediate sheet, a high and usually nonuniform thermal resistance is intro- 
duced and acts in the same manner as insulation between the layers. 

An important conclusion from the above studies is that insulation should 
never be placed between the heat exchanger layers if cooling efficiency is to 
be increased. However, insulation will always increase cooling efficiency if 
appl i ed between the heat exchanger and the hot-gas heat source. Use of the 
upper fin layer as an insulation fin without hydrogen flow will provide equal 
or greater increases in cooling efficiency than folded-in-depth flow for the 
same weight. In fact, a primary reason for rejecting the fo 1 ded- i n-depth flow 
route configuration is that it has two heat exchanger layers, and is thus about 
twice as heavy as single-layer heat exchanger configurations (such as the single 
pass and the folded in width). 

Folded flow wi th multiple injections: Use of folded flow with mu I iple 

injections (fig. 22c) was conceived to reduce cross-section AT' s that occur in 
the folded-in-depth concept and at the same time maintain the beneficial high 
surface temperatures. Some examples of panel systems using this concept are 
shown in fig. 25. Tnis figure indicates the fabrication difficulty, control 
complexity, and, therefore, excessive weight associated with this concept. The 
major limitation on fabrication practicality is the need for injection slots 
with widths of less than 0.004 in. (0.010 cm). 

As noted in fig. 22c, tne excess flow at the right-hand edge of the panel 
may be equal to or greater than zero. The pressure drop is lower, wall and 
hydrogen temperatures are higher, and cooling efficiency is greater for the 
system that has some excess cold flow. 

Flow folded in width: The concept shown in fig. 22d has many counterflow 

streams interspersed across the panel width. A plain-fin or rectangu 1 ar- tube 
configuration is required, since the streams must remain separate between panel 
ends. Use of offset plate fins or pin fins is eliminated, but wavy fins or in- 
serts in the plain passages may be used to increase thermal conductance. This 
concept is called fol ded- i n-wi dth when the inlet and outlet are at the same end 
of the panel (solid lines). It is also possible to have an inlet and outlet at 
both ends of the panel (dashed lines). Both methods of manifolding counterflow 
in adjacent passages develop the same overall performance. The concept with an 
inlet and outlet at both ends of the panel was rejected because of thermal 
stresses caused by a temperature profile similar to that for a single-pass con- 
figuration with a common hot outlet at the panel center and inlets at both ends. 



Initial work with folded-flow panels was based on ac h 'ev.ngwa i tempera 
tures of 2000° R ( I I I I °K) and hydrogen outlet temperatures of 1760 R 1978 kj at 
a heat flux of 10 Btu/sec-ft z ( I 14 kW/m 2 ) , or a hydrogen outlet temperature of 
1 600° R (888°K) at higher heat fluxes. The relative cooling efficiency obtain- 
able with various flow routings was of primary interest. To do this in the 
simplest way, the panel length that could be cooled with a specific flow rate 
and fin geometry was determined. Panels of unequal length resulted. Typical 
temperature profiles and panel length dimensions are shown in fig. 26. With . 
this technique, the flow length at high flux (not shown) was between 4 and 8 in. 

( I 0 and 20 cm) . 

Fig. 26 points out the following characteristics of typical fo 1 ded- i n-w i dth 
performance for panels of unequal length. 

(a) The difference between coolant outlet temperature and maximum cool- 
ant temperature at the hot (folded) end of the panel increases as 
hydrogen Inlet and outlet temperatures increase, but decreases as 
the hydrogen inlet and outlet temperature differential decreases. 

(b) The maximum metal temperature may be reached at the hot end of the 
panel, but further use of the coolant thermal capacity can be 
obtained in a shorter, downstream panel with a smaller hydrogen inlet 
and outlet temperature differential without exceeding the maximum wall 
temperature. Longer panel segments result with higher hydrogen tem- 
perature for the same difference between inlet and outlet temperature. 

A direct comparison of fo 1 ded- i n-wi dth and single-pass concepts was also 
made for panels with a fixed length of 2 ft (0.61 m). This comparison was 
drawn at a recovery temperature of 5000°R (2780°K) and at the hydrogen tempera- 
tures and pressures noted in fig. 27 for high and low fluxes. Results of this 
analysis in terms of cooling efficiency are presented in tabled. Two different 
hydroqen outlet temperatures were used because of differences in the streng o 
superalloys selected for high- and low-flux panels during this part of the study. 
Platn-rectangular fins were used for folded-flow and offset-rectangular f.ns for 
single-pass panels. The Hastelloy X fin and wall (face sheet) thermal conduc- 
tivity was varied as a function of local temperature. A heat exchanger hot-wall 
and cold-wall thickness of 0.010 in. (0.025 cm) was used in all calculations. 

A s i nqle-pass panel was associated in series with a folded-flow panel so 
that identical hydrogen outlet temperatures would prevail for the two concepts 
compared. To obtain maximum coolant outlet temperature, a s . ng 1 e ~ P a s S pane 1 
must be used to end a series of efficient folded-flow panes of equal length. 
Only by use of unequal length or low cooling-efficiency folded-flow panels can 
the Y maximum desired coolant temperature be obtained with folded-flow pane s. 

If a low efficiency folded-flow panel were substituted for the single-pass pa e 
" sel^s wiih the folded panel/the overall cooling efficiency of the system^ 
would increase above the value reported in table 2. The system of two sing e 
pass panels was connected in parallel to allow half the hydrogen to flow through 
each panel. This minimized the pressure drop but did not maximize cooling, 
efficiency. As pointed out for panels with multiple inlets, cooling efficiency 
would havl been increased with part flow through one panel and total flow 
through the other panel. 


For high-flux panels discussed here, all material temperatures were within 
the design limits if the single-pass panel metal temperatures were within these 
limits. The maximum hot-wall temperatures for single-pass and folded-flow 
panels are compared in fig. 28 as functions of flow rate, with fin geometry as 
a parameter, for high flux conditions. The folded-panel system always had lower 
maximum wall temperature than the a 1 1 -s i ng 1 e-pass system. In the system with 
two single-pass panels for heat flux of 500 Btu/sec-ft 2 (5680 kW/m 2 ), maximum 
wall temperature for some fin heights is above the limit for this study: 2000°R 

( M I I °K) . The maximum allowable cold-wall temperature of I600°R ( 888°K) occurred 
at the outlet end of all single-pass panels. 

Fig. 28 also indicates that both single-pass panel and folded-panel cooling 
efficiency improves when fin height is increased. In single-pass panels, cool- 
ing effi ci ency also improves when fin spacing is increased. However, cooling 
efficiency is lowered in folded panels when fin spacing is increased because 
less heat can be transferred from hotter to colder hydrogen. Performance dif- 
ferences between folded and single-pass panels with an increase in fin spacing 
are due to the special function of the hot and cold walls in the folded heat 
exchanger. These 0.0|0-in. (0.025-cm) thick walls act as fins that conduct heat 
from the hotter to the colder hydrogen in adjacent passages. When fin spacing 
is reduced, more heat is transferred from hotter to colder hydrogen, thus increas- 
ing hot-wall temperature and reducing flow rate. 

In fig. 29, inlet pressure is plotted as a function of flow rate and fin 
height for both single-pass and folded-flow systems at high heat flux. The 
folded-flow system always had higher pressure drop and, hence, higher inlet 
pressures than the all-single-pass system. The flow rates quoted for the system 
with a single-pass and a folded-flow panel are one-half the actual flow rate 
per unit width in the panel, and are quoted in this manner for direct compari- 
son with the flow rates quoted for the system with two single-pass panels. The 
actual flow rate in the system with a folded-flow and a single-pass panel is 
twice the quoted flow rate because it cools a 4-ft (|.22-m) total panel length 
compared to the flow rate quoted for the system with two single-pass panels, 
which cools a 2-ft (0.61-m) flow length. Some short fins in the folded-flow 
system have inlet pressure above the design limit in this study of 1000 psia 
(6900 kN/m 2 ). 

The design maximum wall temperature and hydrogen inlet pressure data from 
figs- 28 and 29 were combined in fig. 30 with fin-allowable pressure to show 
that all high flux designs had adequate strength as well as near-minimum weight. 
Minimum weight occurs for all configurations where fin-allowable pressure is 
equal to hydrogen inlet pressure. Minimum weight was a secondary goal in this 
part of the study, since cooling efficiency was known to increase with increas- 
ing fin weight. 

The comparison among all flow routings is summarized in the following re- 
lationships. The pressure drop for single pass is less than for folded-in- 
width, which in turn is less than for fo 1 ded- i n-depth . The cross section tem- 
perature difference for single pass is less than for fol ded- i n-wi d th , which in 


19 



turn is less than for fo 1 ded- i n-depth . These relationships make single-pass 
flow a clearcut choice whenever cooling efficiency is not an important criter- 
ion for panel system design. However, hypersonic vehicle cooling hydrogen flow 
requirements are usually greater than fuel hydrogen flow requirements, so folded- 
in-width flow routing is still of use at low flux and low recovery temperature. 

Temperature profiles: Analysis of flow routing with counterflow in adja- 

cent passages, for both fo 1 ded- i n-depth and folded- i n-width concepts, has indi- 
cated similarities in temperature profile as a function of flow length. The 
hydrogen and maximum metal temperature profiles of single-pass, folded-in-depth, 
and folded-in-width flow routing are compared in fig. 31 at the same hot-gas and 
hydrogen- i n let conditions and for the same panel length. For the folded-flow 
routings, the coolant outlet temperature is lower than both the hotter and 
colder hydrogen stream temperatures along most of the panel. The heat flux 
from the hotter-to-colder hydrogen stream must be greater than the heat flux 
to the panel to produce temperature profiles of this kind. Only with this rela- 
tive flux situation does folded flow increase cooling efficiency. As the ratio 
of aerodynamic heat flux to heat flux from the hotter-to-co I der hydrogen stream 
is increased for the same heat exchanger geometry and hydrogen flow rate, the 
slopes of the hydrogen temperatures are reduced until the average-wall and 
hotter-hydrogen-stream temperatures for folded flow are essentially equal to 
the average-wall and hydrogen temperatures for single-pass flow. This explains 
the fact that cooling efficiency for folded flow is significantly higher than 
for single-pass flow at low aerodynamic heat flux, and negligibly higher at 
high aerodynamic heat flux. 

To maintain an efficient relationship of temperature profiles, acceptable 
cross section AT, and thermal stress at the cold end of folded-flow panels, 
the hydrogen outlet temperature is well below the allowable maximum. Similar 
hydrogen outlet temperatures are used for single-pass flow to provide a realis- 
tic comparison in fig. 31. More than one-half of the available hydrogen thermal 
capacity remains. One or more additional series-connected panels, either single- 
pass or folded, must be used for efficient utilization of the remaining hydrogen 
thermal capaci ty. 

The need to obtain moderate thermal stresses places a severe limitation on 
the cross section temperature difference that is acceptable in all flow routings. 
Fig. 32 shows typical cross section aT profiles that prevail for the temperature 
profiles in fig. 31. At the hot end of the fo 1 ded- i n-depth concept, the cross- 
section aT is equal to the single-pass cross section aT, because there is no 
heat transfer from hotter to colder hydrogen. At the cold end of the folded- 
in-depth concept, the cross section AT is greater than the single-pass AT but 
less than the sum of the single-pass AT plus the local difference between hot- 
ter and colder hydrogen streams. With the fol ded- i n-wi dth flow routing, the 
local maximum surface temperature can be below the local hotter hydrogen stream 
temperature along a substantial length of the panel (fig. 3l). Also, the local 
minimum surface temperature on the adiabatic surface of the heat exchanger can 
be substantially above the local colder hydrogen stream temperature. Therefore, 
the cross section AT for fo 1 ded- i n-wi dth flow routing can be less than 25 per- 
cent of the local difference between hydrogen stream temperatures. 
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High average surface temperature and increased cooling efficiency can be 
obtained in single-pass heat exchangers by use of a thick heat exchanger (tall 
fins) with low thermal conductance and large hydrogen flow area; this obtains 
the large temperature differences primarily responsible for high wall tempera- 
tures. However, unacceptable thermal stresses occur in such a design because 
the low-strength hot end of the panel has almost the same cross section AT as 
the high-strength cold end of the panel. This is in contrast to folded flow, 
which provides a more acceptable small cross section AT at the low-strength 
hot end of the panel, and relatively high cross section AT at the high-strength 
cold end of the panel. 

Insu 1 a t i on . -- Although insulation adds weight and reduces coolant flow 
(increases cooling efficiency) at all recovery temperatures lower than infinity 
the tradeoff between insulation weight and hydrogen flow reduction can be imple 
mented only for a specific vehicle and mission analysis. The efficiency of 
insulation increases as heat flux increases because the insulation thickness 
(weight) required to achieve a given decrease in coolant flow is inversely pro- 
portional to heat flux. Insulation may be powder, batt, foam or solid ceramic, 
a stagnant gas, or, simplest and heaviest, increased heat exchanger face sheet 
th i ckness . 

Fig. 33a depicts a layer of insulating material with a thin refractory 
alloy or superalloy sheet covering, which is in turn held in place by discrete 
attachments. Use of bare ceramic insulation on panels may be limited by the 
need to provide containment, protection from aerodynamic effects, or from 
handling. Oxidation-resistant or coated refractory alloy shields for fibrous, 
powder, or solid insulation are then required. Ox i da t i on- res i s tant coatings 
for refractory alloys are being developed for long term use. Use of oxidation- 
resistant metallic insulation such as Hastelloy X, even though limited to tem- 
peratures around 2500°R (l390°K), was considered because only a thin layer is 
required and it can be included in the brazed panel assembly. 

Fig. 33b shows an overlapping shingle array of metal plates, in which the 
sheets are held in place at their corners. The attachment points consist of 
one fixed support point combined with a set of slotted, oversize holes to pro- 
vide support and yet allow for differential expansion of the shingle with rela- 
tion to the hydrogen-cooled surface. A major defect of this insulation concept 
for hydrogen -coo 1 ed panels is the nonuniformity of thermal resistance between 
the shingles and the heat exchanger face sheet. This nonuniform resistance 
can cause uneven shingle temperatures and oxidation rates or local hot spots 
in the heat exchanger, which are detrimental to cycle life through increased 
thermal stresses. The surface roughness inherent in the overlapping shingle 
concept will be a disadvantage in some applications where smooth aerodynamic 
surfaces are required. 

A metallic sheet that is held away from the hydrogen-cooled surface by pin 
fins is shown in fig. 33c. The pin fins hold the metallic sheet and provide a 
uniform heat path with a good thermal resistance to weight ratio. The face 
sheet can be made up from many small elements where the differential expansion 
between them and the hydrogen-cooled surface can be absorbed by lateral bending 
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of the pins in any direction. The use of the same cross section area in plate 
fins of either the plain or offset rectangular type can provide the same insu- 
lation effects for equal weight, but will have bending flexibility in one 
direction only: perpendicular to the corrugations. 


A recovery temperature of 5000°R (2780°K) was selected as typical for many 
hydrogen-cooled panel applications to illustrate the effects of meta 1 1 i c i nsu- 
lation. Typical insulation performance is listed in table 3, as obtained by 
the method described in Appendix A. The tabulated hot-gas heat transfer coef- 
ficients were used for both insulated and uninsulated surfaces because the hot- 
gas heat transfer coefficient is a weak function of wall temperature.. The heat 
fluxes resulting at the hot end of the panel with and without insulation are 
for the tabulated local wall temperatures. A realistic maximum temperature for 
nonstructural superalloys of 2500°R (l390°K) permitted a panel hot-end heat 
flux reduction of 20 percent. In this case, the hydrogen flow rate and average 
heat flux were reduced more than the local heat flux at the panel hydrogen out- 
let end, because superalloy thermal conductivity increases with increasing 
temperature. The flow rates in table 3 apply for the integrated heat flux on 

the entire panel length. 

A Haste 1 loy X hot wall with a thickness greater than 0.010 in. (0.025 cm) 
can provide the required insulation but is heavy. Insulation weight can be 
reduced using a fin layer and a thin face sheet instead of solid metal. The 
thin insulation closure sheet provides a smooth aerodynamic surface, and con- 
ductivity of the stagnant gas in the fins will be negligible compared with 
conductivity of the metal fins. For the typical sheet and fin insu^tion 
geometry analyzed under typical conditions; less than 0.7 Ib/ft (3. g m 
insulation can reduce hydrogen flow rate more than 20 percent. 

Typical fins, such as 20/in. (7.9/cm) with a thickness of 0.006 m. (0.015 
cm) or 40/in. ( 15. 8/cm) with a thickness of 0.003 in. (0.008 cm), have a soli- 
dity of 0.12 and an effective thermal conductivity about 8.3 percent that of 
solid metal. Use of equation 10 in ref. 5 indicates that air conduction and 
radiation will increase the effective fin conductivity noted by less than 3 
percent and I percent, respectively, at the hydrogen outlet end of the panel. 
Therefore, the average effective thermal conductivity of the fins described 
above is 8.6 percent of solid metal conductivity, and is practically equal 
the effective thermal conductivity based on metal conduction alone. 


Insulation fin temperature differences are less than the overall 
insulation temperature differences because the insulation hot wall has a 
significant AT at high heat flux. The insulation fin height provides for 
a length of metal between fin roots equal to that required for the insulation 
fin AT, since fin thickness was added to wall thickness to obtain wall AT. 
Insulation fin heights are shorter than those normally used in fluid passages, 
but the heights cited have been fabricated. Weight and thermal resistance of 
braze material have not been included in the tabulated values. Since insula- 
tion thickness is based on allowable maximum surface temperature at the hot 
end of the panel, the tabulated weights per unit area are independent of 
overall panel dimensions and depend only on the density of the insulation fin 
material. A density of 0.3 lb/in. 3 (830 kg/m 3 ) was used. 
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Manifold Performance 


The objective of the manifold study was to provide estimates of pressure 
loss and flow nonuniformity resulting from various manifold concepts for use 
in panel concept evaluation and tradeoff studies. A low manifold pressure 
drop is desirable since the manifold contributes to the overall pressure drop 
of the cooled panel. However, pressure drops greater than the minimum attain- 
able may be required to provide proper flow distribution with low manifold 
weight. For a uniform heat flux (which was generally assumed in this study), 
a uniform flow is desirable to prevent nonuniform temperatures. The effect of a 
typical type of flow nonuniformity Is shown in fig. 34, where 10 percent more 
flow occurs in the center of the panel than at the panel edges. The temperature 
is higher along the edges of the panel than in the center. This nonuniform tem- 
perature results in (l) thermal stresses than can easily attain large values, 
which must be allowed for in design, and (2) overheating of the starved areas of 
the panel. Since maximum temperature capabilities of materials are limited, it 
becomes necessary to increase the total coolant flow by 5 percent, for this 
example, to prevent overheating (ref. l). However, thermal stresses attributable 
to temperature nonuniformity In the panel width are not eliminated by simple 
increases in coolant flow. With uniform heat flux, increases in thermal stress 
and coolant flow rate are avoided only with uniform flow. Thus, correct mani- 
fold design provides a match between cooling requirements and coolant flow dis- 
tribution across panel width to minimize thermal stresses and coolant flow rate 
requi rements. 

Manifold concepts suitable for the special requirements of regenerat i ve 1 y 
cooled panels are presented in ref. I. Three basic concepts, shown in fig. 35, 
were selected for use in this program. The upper corrugation or fin is the 
heat exchanger core fin in all of the concepts depicted. The two types of 
flat manifold provide pressure containment by fins similar to those used in 
the heat exchanger core. Cylindrical manifolds of the type shown can be used 
to obtain low pressure drop and uniform distribution, but result in relatively 
large unsupported spans that must accommodate high coolant pressures. 

The selection of these three basic manifold concepts was based on the need 
for assembly of adjacent panels with the use of edge seals to prevent hot gas 
flow into regions in back of the panel. Additionally, the concepts provide for 
coolant flow and high thermal conductance in all parts of the panel surface 
exposed to hot gas. 

Basic assumptions and constraints . — The flat rectangular manifold 
(fig. 35a) was selected as a reference configuration, with uniform flow (zero 
maldistribution) as the design objective. The rectangular manifolding was 
chosen to take advantage of previous design layout and analysis work accom- 
plished for this configuration. It was assumed that the trends established 
from analyzing the rectangular manifold would be representative of the tapered 
manifold also. Consequently, pressure drop calculations were performed only 
for the rectangular manifolds. The cylindrical manifold was not considered 
beyond the initial stage of the study because the concept has a definite weight 
disadvantage for containing coolant pressure. 
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It is possible for uniform flow to occur through a cooled panel, indepen- 
dent of the properties of a manifold, if the panel pressure drop is very large 
compared to the manifold pressure drop. When this occurs, any nonuniformity 
in the pressure drop through the manifold is negligible compared to the uniform 
resistance across the width of the panel heat exchanger core. For this mani- 
fold study, however, it was assumed that heat exchanger pressure drop would be 
low enough that a variation in pressure drop through the manifold would have an 
effect on flow distribution. Thus, if uniform flow is achieved in the manifold, 
the ratio of heat exchanger panel pressure drop to overall pressure drop can 
have any value, and uniform flow will still occur. 


Uniform flow in the manifold dictates that the resistance between the 
manifold inlet port and the heat exchanger panel be equal for all paths. 
Therefore, all paths except the maximum length paths that feed the edges o 
the panel require some source of pressure drop in addition to the pressure 
containment fins. Two methods for obtaining equal pressure drop for all flow 
paths are: (l) variation of the resistance in the vertical fin (the fin 
between the upper edge of the port and the heat exchanger core fin) an 
(2) incorporation of an orifice plate with large free-flow area in series 
with the longest, most highly restrictive paths at the edges, and small flow 
areas in series with the shortest, least restrictive paths adjacent to the 
ports. This study considers the second method. 


The large number of dimensional variables involved in manifold design 
require selection of additional guidelines and assumptions for size and con- 
struction in order to facilitate calculation of pressure drop and weight. 
Accordingly, the pressure drop analyses were based on the following geometric 
features shown conceptually in fig. 35a: 


(a) The configuration is flat rectangular. 

(b) One manifold port feeds the full heat exchanger width. 

( c) The manifold port and associated piping diameters are sized to 
provide a ratio of core free-flow area to port free-flow area 
of I. Specification of this ratio provides a convenient basis 
for sizing the piping relative to the heat exchanger. 

( d) The overall vertical fin flow length of 1.5 in. (3.81 cm) is divided 
into two parts: 0.5 in. (1.25 cm) from the top edge of the port to 
the bottom surface of the panel, and I in. (2.54 cm) parallel to the 
panel surface and underneath it. This I -in. (2.54-cm) length con- 
tains three 90-deg bends and accommodates the seals, preventing hot 
gas flow to the region behind the panels. These di mens i ons are 
typical values and were selected with the aid of layout design work 
presented in ref. I. 

( e) Pressure containment in the manifolds results from plain rectangular 
fins, 10 per in. (3.94 per cm) of 4-mil (0.010-cm) Hastelloy X and 
manifold face sheets of 0.010-in. (0.025-cm) Hastelloy X. These 
selections are adequate for pressure containment of 500 psi 

(3450 kN/m 2 ) at I600°R (888°K) and represent minimum gauge limits. 
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( f) Manifold fin heights of 0.250 and 0.050 in. (0.634 and 0.127 cm) are 
used to provide an adequate design range for both high and low flux 
pane 1 app 1 i cations. 

(g) Inlet and outlet manifolds are identical in size, although typical 
design point densities require that the design be based on the outlet 
manifold, which contributes 97 percent of the total manifold pressure 
drop. When inlet and outlet manifolds are the same size, the weight 
penalty is small, and fabrication is simplified. Also, off-design 
operation benefits from equal manifold sizes. 

(h) The port-shoulder-radius to port-diameter ratio is 0.1 at the transi- 
tion point between the inlet or outlet pipe and manifold. This value 
was selected to provide a calculation basis, although its effects are 
only approximated in the calculation procedure. 

Analytical results . --Using the assumptions and constraints just discussed 
and procedures described more fully in Appendix C, pressure drop and weight were 
calculated for specific manifold configurations. Manifold performance was 
related by a flow rate parameter that combines heat flux and panel dimensions 
and is proportional to hydrogen flow rate, specific heat, and temperature change. 
This parameter, (q/A)(^/w), permits flow rate to be indicated without a separate 
curve for each flow rate in each panel width. Fig. 36 shows the combined inlet 
and outlet manifold pressure drop as a function of the flow rate parameter for 
three different manifold widths and port diameters. The results are for an inlet 
coolant temperature of I00°R (55.5°K) and an outlet temperature of I600°R 
(888°K). Other inlet and outlet hydrogen temperatures will cause a negligible 
difference in oAP. This is because manifold friction pressure drop is a weak 
function of viscosity, and friction does not exceed one-half the manifold pres- 
sure drop. The outlet pressure used to calculate the pressure drops shown in 
fig. 36 was 250 psi (1560 kN/m 2 ). Consequently, the curves in fig. 36 indicate 
pressure drops slightly higher than these that occur with high flux panels and 
slightly lower than those that occur with low flux panels because the inlet 
pressure used in developing the data was 500 psi (3120 kN/m 2 ). The inlet mani- 
fold pressure drop is about 3 percent of the combined inlet and outlet manifold 
pressure drop at the pressure and temperature conditions that form the basis for 
fig. 36. Inlet manifold pressure drop, although linearly proportional to hydro- 
gen inlet pressure, Is considered constant because it is such a small fraction 
of the combined inlet and outlet manifold pressure drop. 

The limiting value of the flow rate parameter at any particular manifold 
width Is noted when the port pressure drop is equal to the overall manifold 
pressure drop allotment. This limit requires a zero pressure drop, infinite 
fin-height manifold. Work in support of the studies in ref. I, for example, 
used a pressure drop allotment of 45 psi (310 kN/m 2 ) for the manifold and 
assumed that 40 percent of the manifold pressure drop occurred in the port. 

On this basis, the manifold pressure drop allotment is equaled by the port 
pressure drop when the combined inlet and outlet manifold pressure drop is 
112 psi (770 kN/m 2 ). Fig. 36 shows that the limiting value of the flow rate 
parameter is 600 Btu/sec-ft 2 (6810 kW/m 2 ) for a manifold width of 2 ft. 
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For flow rate parameters greater than the above limit, two options are 
available. One is to increase the port diameter and the other is to reduce 
manifold width by increasing the number of ports in the panel width of 2 ft 
(.61 m). The option of more than one port per panel width was exercised. 

Table 4 lists the selected number of manifold width segments and ports, mani- 
fold fin height, pressure drop, and port diameter. Manifold widths of 8 and 
12 in. (20.3 and 30.5 cm) were used for the 500 Btu/sec-ft 2 (5680 kW/m 2 ) heat 
flux with a 3-ft (.9l4-m) flow length to illustrate the change in manifold fin 
height that occurs. The 8-in. (20.3-cm) segment width, which yields a shorter 
manifold fin height, was used for subsequent analysis. 

Calculated manifold weight as a function of manifold width is shown in 
fig. 37 for tapered and rectangular manifolds. The port diameter is plotted 
because of the fixed ratio of port diameter to manifold width. Tapered mani- 
folds are lighter than rectangular manifolds because the 0.5-in. (1.27-cm) 
vertical fin length between the top of the port and the bottom surface of the 
panel was tapered. The possible weight reduction afforded by use of tapered 
manifolds is most important at low heat-flux, low normal-pressure load design 
points, where the manifolds are a larger fraction of the total panel weight 
(ref. l). Fig. 38 combines results from figs. 36 and 37 to show the trade 
between: (l) installation and ducting complexity, and (2) manifold weight and 

pressure drop. As the manifold width fed by single duct is increased to 
simplify installation, the manifold pressure drop and/or weight increases. For 
the geometry assumed, as manifold width increases, pressure drop increases much 
faster than manifold weight. Other assumptions for geometry (larger fin 
height and port diameter) can provide wider manifolds with no increase in 
pressure drop but at a higher rate of weight increase. 

The manifold pressure drops shown in fig. 36 are based on uniform flow 
distribution across the width of the manifold. As already discussed, this is 
a convenient and practical approach to the manifold design problem. Fig. 39, 
however, shows the type of flow distribution obtainable in a typical manifold 
when no special provisions are made to obtain uniform flow. Table 5 shows the 
manifold geometry and operating conditions used in calculating the flow distri- 
bution of fig. 39. The pressure drops associated with the three widths of 
manifolds are shown in fig. 40. The maxi mum-to-mi ni mum flow rate ratios and 
the ratio of core pressure drop to overall pressure drop are also noted. As 
shown in the experimental results, the core-to-overal 1 pressure drop ratio is 
the most useful parameter for providing an indication of flow distribution in 
terms of maxi mum-to-mi ni mum flow ratio. Only with the relatively high pressure 
drop ratios noted in fig. 40 can the relatively small flow nonuniformities 
indicated in figs. 39 and 40 be achieved* 

Exper i mental resul ts . — A completely analytical treatment of the manifold 
flow distribution problems has not proved successful in the past, even for 
relatively simple configurations. The purpose of tests performed during this 
program was to evaluate the performance of manifold configurations considered 
practical for application to regene rat i vel y cooled panels. A widely applicable 
configuration is the flat manifold, rectangular or tapered, attached to the 
regene rat i ve 1 y cooled surface at right angles (figs. 35a and 35b). This was 
the reference configuration selected for evaluation. Table 6 lists the test 
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specimen geometries evaluated with the referenced configuration. Dimensions 
shown were obtained by measuring the parts. The test apparatus and procedures 
are described in Appendix D. 


The maximum-to-mi nimum flow ratio through the core (the regenerati ve 1 y 
cooled surface) is the characteristic flow uniformity parameter. Preliminary 
analysis indicated that the geometric manifold parameters of interest were the 
'? 0r port dlam f ter > plate spacing (manifold fin height), manifold 

width parallel to panel width, and port radi us-to-di ameter ratio. Figs. 41 and 
42 summarize the experimental investigation of the effect of these geometric 
parameters on flow uni formi ty and on pressure drop. Appendix D gives details 
o t e ata analysis as well as detailed performance data for the various test 
specimens. In general, the data reveal that increased flow uniformity results 
from manifold geometries with reduced pressure drop (resistance), because 
resistance is unequal in the various manifold flow paths. 

In addition to pressure drop ratio, geometric similarity also provides a 
good first-order indication of flow distribution uniformity. Specimen 12 is 
twice the size of Specimen 2 in manifold width, port diameter, fin height, fin 
spacing, and fin thickness. Rather close agreement in flow rate ratio is 
obtained with W max /W m ; n of 1.88 for Specimen 2 and 1.69 for Specimen 12. 

Fig. 41 shows that flow is more uniform as port diameter is increased, as 
mani fold wi dthi s reduced, as port radi us-to-di ameter ratio is increased, and 
as fin height is increased. Although changing port r/D provides as much change 
in flow uni formi ty as port diameter, the reason for this is not primarily a 
change in port loss coefficient but a change in horizontal fin friction and 
flow-turning pressure drop between the horizontal fin and vertical fin. With 
large port diameter or large r/D, the width of the horizontal fin adjacent to 
the port is approximately the same. Specimens 4 and 7 have practically equal 
flow nonuniformity (fig. 41) and horizontal fin flow width, but have markedly 
different port diameters (table 6). Even for the relatively narrow spans ' 
evaluated as part of the experimental program, flow distribution ratios (maxi- 
mum ow to minimum flow at any station) of up to 5.75 were experienced 

[ f '£\ i\{’ 0n ^ reason for this is that the configurations were selected so as 
to highlight the problem rather than to minimize it. For example, the manifold 
plate spacing (fin height) for most of the specimens was set at 0.052 in. 

(0.132 cm); this resulted in relatively severe pressure drops and large mal- 
distribution. In general, the weight penalty associated with greater plate 
spacings is small, so that use of such plate spacings appears quite feasible. 

Use of isothermal air was a second characteristic of the experimental evalua- 
tion that exaggerated flow nonuniformities. Both inlet and outlet manifolds 
contributed significantly to flow nonuniformity. For the boundary conditions 
used in the heat transfer analysis work of T r = I00°R (55.5°K) and T 

1600 R (888°K), flow distribution is largely controlled by the lower density 
hydrogen in the outlet manifold. 


The lowest test value for maxi mum-to-mi n imum flow ratio of 1.32 will 
probably not be acceptable for any operational application because of hydrogen 
flow rate and stress increases. Extrapolation of the curves to higher ratios 
of core to overall AP indicates that maxi mum-to-mi nimum flow ratios of l.l or 
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less require a core to overall-AP ratio greater than 0.5. Uniform flow wi 1 1 
occur at a AP ratio of 1.0. In general, both analysis and test results indi- 
cate that manifolds for use in regenerati vel y cooled panels, even where the 
ratio is above 0.5, will require means for equalizing re si stance in all 
manifold flow paths. 

The experimental data can be used to make first estimates of flow distri- 
bution for configurations other than those tested. Appropriate correcting 
devices can then be designed for experimental evaluations. Any correcting 
device for configurations of the complexity required for panel manifolding 
will require iteration of the design based on test results. 


CONCLUDING REMARKS 


The heat transfer and fluid flow aspects of hydrogen-cooled flat panels 
have been evaluated. Pressure drop and thermal conductance data for a wide 
variety of coolant passage configurations were obtained. Application of data 
to integrated panel design was undertaken and is illustrated for typical cases 
in the present report. 

Selection of a heat transfer design for a given application generally 
requires consideration of other factors. These factors include weight, thermal 
stress, structural life, system pressure levels, coolant consumption, manu- 
facturing methods, and installation. No single performance figure of merit 
has been found that permits selection of coolant passage geometry, even if 
heat flux, thermal stress, coolant consumption, and coolant pressure drop 
are retained as the only parameters. Recourse to basic thermal conductance 
and pressure drop data, on the other hand, can quickly serve to narrow the 
selection to a limited and manageable number of candidates. 

Various schemes for conserving coolant were investigated. Insulation is 
a logical candidate at all times for reducing heat load to the coolant. The 
benefits of flow folding by various techniques were also evaluated. As with 
insulation, flow folding aims at raising the average temperature of the heated 
surface and, hence, reducing heat load. At low heat fluxes, and at recovery 
temperatures below 7000°R (3890°K), cooling efficiencies can be increased by 
5 to 10 percent. As heat flux and recovery temperature increase, the benefits 
are reduced. The complexity associated with any of the flow folding schemes 
can, therefore, be justified only where large areas operating at low heat-flux 
levels constitute a substantial part of a total flight system heat load. Pre- 
liminary analysis of the problem, prior to availability of the detailed analy- 
sis program, indicated much more substantial benefits from flow folding. The 
final results obtained were, therefore, somewhat disappointing in view of 
expected benefits. 

The co nee rn in typ ical panel applications withe oo 1 an t conse rvat i on , 
thermal stress, pressure drop, envelope, and weight requires consideration 
of not only the panel proper, but of the manifolds. Indeed, the benefits 
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obtained by careful optimization of the coolant passage geometry and by the 
use of insulation and flow folding can be completely obviated by relatively 
small amounts of flow maldistribution in the manifolds. 

The manifold studies presented in this report provide pressure loss and 
flow nonuniformity data resulting from use of various mani fold concepts. 
Extension of this data to configurations other than those specifically evalu- 
ated is feasible, but will generally require experimental verification and 
iteration of the design. As a starting point for design, core ( panel )-to- 
overall pressure drop ratios of 0.5 or higher appear necessary for satis 
factory manifold performance, with maximum-to-mi ni mum flow ratios of l.l or 
less across the width of the manifold. Even at these pressure drop ratios, 
efficient manifold design, i.e., light weight and low volume, may require use 
of special inserts or orificing to equalize pressure drops and flow rates 
across the panel. 

All performance data presented in this report are based on the use of 
fluid properties evaluated at coolant bulk temperatures. The selection of 
reference temperature and basic performance correlation can significantly 
influence predictions relating to panel wall temperature difference. This 
is particularly true where wall-to-bulk temperature ratios are large and 
where hydrogen temperatures are below 90°R (50°K)._ The type and form of 
correlation used in design must, therefore, be an important consideration 
in regene rat i ve 1 y cooled panel applications. 


29 


always used. For curves with heat flux as 
in hydrogen enthalpy between I00°R (55.5 K) 
(3450 kN/m 2 ), i.e v W/w= 0.0001477 ( q/A)A , 
kg/sec-m) . 


absci ssa, W/w was based on the change 
and 2000° R (MM°K) at 500 psia 
1 b/sec-f t (W/w = 0.00022 (q/A)^, 


Heat exchanger thermal conductance, heat flux, and fin AT are related by 
eq. 3. 


T W - T F = aT fin 


T _ T . 9 Lh. 

W C hAj / vd 


(3) 


A 10-percent margin for manufacturing and curve-reading tolerances is included 
■„ alfestiltes of T,„hA ,/wt and T w - T f . That is, values of T.hA./wA presumed 

in this report are 0.9 of the calculated values and values of T w - T p given are 
l.ll of the calculated values. 

When the hot wall was evaluated as insulation, k. was evaluated at 

(T + T )/2, where T = 2500°R ( I390°K) and T y = 2000°R ( I I I 1 0 K) ; i.e., panel 

hotfend k , = k. ,. =16.2 Btu/hr-°R-ft (28 W/m-°K) at 2250°R (|250°K) 

nut enu i'.. p | f | n 

Use of a fin layer brazed under a thin cover sheet was also analyzed as 
insulation. The above constants, the values for panel hot-end heat ux rom 
table 3, and an insulation AT of 500°R (278°K) were substituted in eq. * to 

solve for h. fJn . The insulation thickness and weight were then calculated from 

eqs. 5 and i.^Eq. 4 was developed for the assumption of an effective insulation 
wall thickness of t. f + t. fJn and an effective insulation fin height of 


h. 
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Single-Pass Heat Exchangers 

All single-pass heat exchanger thermal conductances and pressure drops at 
infinite recovery temperature were calculated by a digital computer program 
using the method described below. 
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Except where noted, the heat exchanger thermal conductance data are for 
local conditions of hydrogen temperature and pressure at the hydrogen outlet 
(hot end) of the heat exchanger. The calculation procedure shown below is also 
for local conditions, since overall performance for panel heat exchangers is 
usually obtained by integration of a series of incremental calculations, each 
made at local conditions. Specific assumptions are listed relative to appro- 
priate parts of the analysis. 

The cal cul at i on procedure described here is based on pipe flow correlations 
and applies to single-phase flow only. Conclusion 2 in ref. 8 indicates the 
general correctness of using pipe flow equations with all fluid properties eval- 
uated atji 1m temperature when, as in this study, bulk temperatures are above 
90 R (50 K) . Ref. 2 suggests the use of bulk temperatures in frequently inter- 
rupted boundary layer heat exchanger geometries, although test data at large 
wall-to-bulk temperature ratios are lacking. Film temperature in this discussion 
is equal to the arithmetic average of wall and bulk temperatures. The general 
curves of thermal conductance and fin AT were developed with fluid properties 
from fig. 2 evaluated at 2000°R ( I I II 0 K ) , a value that may be considered bulk 
or film at the user's discretion. If the temperature is considered to be film 
temperature, the flow rate parameter, W/w, must be multiplied by the ratio of 
bulk-to-film temperature before reading thermal conductance. 

Bulk temperature was used to obtain all of the performance curves with 
cross section AT or wall temperature. It is convenient to use the local bulk 
temperature because iteration is avoided. Iteration is required with use of 
film temperature or a wall-to-bulk temperature ratio. Use of bulk temperature 
and wall-to-bulk temperature ratio results in nearly the same hot-wall tempera- 
tures at 1600 R ( 888°«) hydrogen temperature. Wall temperature and fin AT are 
far higher when wall-to-bulk temperature ratio is used at a I00°R (55.5°K) 
hydrogen temperature. Use of wall-to-bulk temperature ratio gives essentially 
the same results as use of film temperature when all fluid properties are evalu- 
ated at film temperature, and the heat transfer coefficient is multiplied by the 
bulk-to-film temperature ratio to the 0.8 power. 


The parameter j is defined by eq. 7 and read from curves as a function of 
Reynolds number (eq. 8). 


j 


St (Pr) 2 / 3 


h A f (Pr) 2 / 3 


W C 

P 


Re 


4 r.W 
h 

M-A f 


(7) 

( 8 ) 


These j vs Re curves were developed from data on the same or geometrically 
similar fluid passage geometries reported in ref. 2. The specific figure from 
ref. 2 used for f and j data is noted as the last term in the various fin 
geometry designations. The fin geometry designations are defined in the 
Symbols and Parameters table. 
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Specific test data are available for the 16 fin/in. (6.3 fin/cm) and the 
20(7. 9)R - 0.100(0.254) - 0.125(0.318) - 0.004(0.010) fin geometries. Fig. 

10-61 from ref. 2 was used for all other rectangular offset f i n geometri es. To 
make use of a single performance curve suitable for a broad variety of fin 
geometries, the ratio of offset length, L q , to hydraulic radius, r h , was held 

constant at 8. The strongest control of the level of j and f for a particular 
class of fin geometry is the ratio of offset length to hydraulic radius. The 
constant 8 is empirical; it is based on test data in ref. 2 for geomet ri es wi th 
offset length of 0.125 in. (0.318 cm) or greater, 20 fins/in. (7.9 fins/cm) or 
less, height of 0.100 in. (0.254 cm) or greater, and thickness of 0.004. in. 
(0.101 cm) or greater. Fig. 43 shows relationships between spacing, height, 
thickness, offset length, and hydraulic radius for rectangular fins of interest 
in this study, for which performance was calculated by fig. 10-61 of ref. 2. 

Eq. 9 was used throughout this study to calculate heat transfer coeffi- 
cients. Eqs. 10 and II are widely used correlations that evaluate fluid proper- 
ties at other than bulk temperature. 


bulk h = 


j W c 


A f ( Pr) 2/i 


(9) 


where all fluid properties are evaluated at local bulk temperature. 


jf W C D f ^ T bu 1 k /T f ^ 
fi1m h = £ (Pr f ) 4/i 


0 . 8 


(io) 


where all fluid properties are evaluated at the arithmetic average of wall and 
local bulk temperature. 


wal 1 -to-bul k h = 



(ID 


where all fluid properties are evaluated at local bulk temperature. For 
Reynolds numbers less than 5000, n is 0.0; for Reynolds numbers greater than 
5000, n is 0.55. 

Overall heat transfer surface effectiveness is calculated by eq. 12 for 
fins with uniform thickness parallel to flow. 
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where for rectangular corrugated fin geometry 

A fin =«( 2 (N)(h fIri -t f . n ) + (l ' 'fin N >) 
a t = A fin + w«(l - t fin N) 

^fin = ^ Btu/hr-°R-ft (17.3 W/°K-m) unless otherwise noted 


1 r . “ hr- + 0 

fin fin 


• 5 (i!r- 3 t fin) 


A f = w(l -tj.. N)(h.. - t _ . ) 

f fin fin fin ; 


h A f 

r h = — T (= D/4 for round tubes) 


( 13 ) 

(14) 

( 15 ) 

(16) 

(17) 


Friction pressure drop was calculated by eq. 18. The Fanning friction 
factor, f, was read from curves as a function of Reynolds number (eq. 8). 
These f vs Re curves were developed from data on the same or geometrically 
similar fluid passage geometries. 


AP 


f 


f l W 2 

V 2g c p a 


( 18 ) 


Hydrogen properties (fig. 2) at I050°R (584°K) and 500 psia (3450 kN/m 2 ) were 
used to estimate ctAP/£> oAP, and AP. A 10-percent margin for manufacturing and 
curve-reading tolerances is included in the stated values of oAP/j£, resulting 
in l.l times the calculated values. 


Pressure drop due to change in momentum, calculated by eq. 19, was not 
included in any of the aAP curves, but was included In the curves for p . A 

constant ratio of AP m /AP^ = 0.06 was used. 


AP 


m 



(19) 


Flow Routing Analysis Methods 

All flow routings herein discussed, including the single-pass configuration, 
were analyzed by means of digital computer programs. Fig. 22 shows that simul- 
taneous heat transfer can occur between the hot gas and all coolant streams only 
in the fo 1 ded- i n-wi dth concept. The other folded-flow configurations have heat 
transfer between the hot gas and only one coolant layer. Separate analysis pro- 
grams were used for the two general cases. 
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Both programs permit only two of the fluids to have temperature changes. 

The other one or two fluids must have constant temperature, which implies 
infinite capacity rate. This occurs for some two-phase heat transfer processes, 
and is approached for aerodynamic heating at hypersonic flight conditions. The 
heat transfer coefficients for the infinite capacity rate fluids are input as 
constants. Passage geometry, flow rate, temperature, and pressure are assumed 
to be uniform in the flow width. Incremental analysis is used, with the maximum 
length and incremental length selected as inputs. 

A three-fluid heat exchanger analysis program was used for the folded-in- 
depth concept, and is applicable to plate and rectangular-fin heat exchanger 
geometries, including those having no fins in either one or both of the edge 
fluids. The finite capacity rate fluids (central and one edge) may be any 
single-phase liquid or gas. One edge fluid must be of infinite capacity rate 
because it is analyzed as rejecting or absorbing heat without temperature change . 
A range of fluid flow rates and boundary temperatures was used with the fin 
geometries considered to find the combination that gives the desired heat 
exchanger length. 

The heat exchanger geometry inputs are: 

Metal thermal conductivity 

Fin spacing, fin height, and fin thickness for each of three fluids 

Initial number, final number, and increment in number of fin layers 
(sandwiches) for each of three fluids in each three-fluid cycle 

Initial width, final width, and increment in width of heat exchanger 

Initial number of cycles, final number of cycles, and increment in 
number of three-fluid cycles 

Initial length, total heat exchanger length, and length increment 
between initial and total length at which calculations are to be 

made 

Colburn modulus, j, and Fanning friction factor divided by Colburn 
modulus, f/j, as functions of the Reynolds number for each of two 
fluids undergoing temperature change 

The operating condition inputs are: 

Flow rate for both fluids undergoing temperature change 

Constant temperature for infinite-capacity-rate fluid and temperature 
at some point in the heat exchanger for the two fluids undergoing 
temperature change 

Heat transfer coefficient for infinite-capacity-rate fluid 
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Specific heat, density, dynamic viscosity, and (Prandtl 
for both fluids undergoing temperature change 


number) 


2/3 


/C 


P 


The outputs are: 


Input values of fluid flow rates and passage geometries 

Calculated values of Reynolds number and heat transfer coefficient 
for the two fluids undergoing temperature change 

A tabulation of fluid and wall temperatures, heat fluxes, temperature 
differences, and density-adjusted pressure drops as a function of 
position relative to one end of the heat exchanger for both a three- 
fluid heat exchanger and a two-fluid heat exchanger 

The two-fluid heat exchanger is similar to the three-fluid heat exchanger 
except that in the former, the outer fin layer with a finite-capacity rate is 
deleted. This feature permits direct comparison for identical geometry and 
fluid inlet conditions of both a fol ded- i n-depth and a single-pass heat 
exchanger in one computer output. 


The analysis procedure for each length increment requires calculation of: 

(a) The flow passage geometry for each fluid in terms of free-flow area, 
hydraulic radius, fin heat transfer area, and plate heat transfer 
area from input dimensions 

(b) Heat transfer coefficient, thermal conductance, and aAP for each of 
the two fluids undergoing temperature change from pipe-flow equations 
and j and f/j tables 

( c) Wall temperatures for the center fluid passage from a heat balance on 
these walls and the one-dimensional conduction equation applied to the 
fins connecting these walls 

( d) Temperatures of the finite-capacity-rate fluids at the end of the 
length increment (used as inlet to the next increment) by two simul- 
taneous differential equations based on a heat balance of the finite- 
capac i ty- rate fluids and the wall temperatures from procedure ( c) 

(e) Net heat fluxes and the other wall temperatures by use of the values 
obtained in steps (b) and ( c) above 


A four-fluid heat exchanger analysis program was used for the folded-in- 
width concept. A cross section and one flow configuration analyzed by this 
program are depicted in fig. 44. The four fluids are identified as G, C, H, and 
A. Flow directions of fluids G and A are unimportant because both fluids are 
i nf i ni te-capaci ty-rate ; that is, each fluid has a constant temperature, whether 
heat is being added or removed. The heat transfer coefficients for fluids G 
and A are assumed to be constant and are inputs to the program. Fluids C and H 
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are in passages between fluids G and A, and may be any single-phase liquid or 
gas. They may flow parallel or counter to each other. 

For fluids C and H, the input includes: 

Flow passage geometry (cross section dimensions and metal thicknesses) 

Tables of Fanning friction factor, f, vs Reynolds number, Re 

Tables of Colburn modulus, j, vs Reynolds number. Re 

Tables of fluid properties (specific heat, dynamic viscosity, and 
Prandtl number, all as functions of temperature) 

Flow rates, WC and WH 

Inlet temperature for fluid C, TC, and outlet temperature for fluid 
H, TH 

Flow geometry (parallel or counterflow, folded or not folded, total 
length, and analysis increment size) 

A table of metal thermal conductivity as a function of temperature 
For fluids G and A, the input includes: 

Temperature and heat transfer coefficient 
The output includes: 

All inputs of importance 
At each station of interest: 

20 metal and 2 fluid temperatures (locations indicated in fig. 44) 

Heat flux and heat transfer rate to or from surfaces facing 
f 1 ui ds G and A 

ctAP for fluids C and H 

Reynolds number, heat transfer coefficient, and heat transfer 
rate to or from fluids C and H 

Although the input includes dimensions and flow rates for one complete 
passage of fluids C and H, the temperature distributions calculated are for 
the shaded area between the lines of symmetry in fig. 44. The analysis proce- 
dure for each increment involves calculation of: 

(a) The free-flow areas, heat transfer areas, and hydraulic radii of the 
passages for fluids C and H from input dimensions. 
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(b) The heat transfer coefficient, thermal conductance, and aAP for 

fluids C and H from pipe-flow equations, using j and f tables. 

(c) A heat balance at nodes I and 2 to obtain Tl and T2 in terms of the 

local fluid temperatures, heat transfer coefficients, and geometry. 

As a basis for thi s cal cul ation, a one-dimensional fin temperature 
distribution analysis is applied to the partition fin in terms of 
Tl and T2, TC, TH, and the heat transfer coefficients for fluids C 
and H. Also, a one-dimensional fin temperature distribution analysis 
is applied to passage walls acting as fins for the assumptions that 
all values of Tl are equal and that all values of T2 are equal. Heat 
transfer coefficients and temperatures for fluids A and G are' used in 
this part of the analysis. 

( d) Temperatures TC and TH at the end of the increment (the inlet to the 

next increment) by two simultaneous differential equations based on 

a heat balance on fluids C and H. 

(e) The five temperatures along the centerline of each passage wall from 
temperatures Tl and T2, and the one-dimensional fin temperature dis- 
tribution analysis mentioned in ( c) above. 

( f) The 10 external surface temperatures and eight internal surface tem- 
peratures for the passage walls between the f i n i te-capac i ty- rate 
fluids C and H and the infinite-capacity-rate fluids G and A, from a 
one-dimensional conduction analysis using the local fluid temperatures 
and heat transfer coefficients as well as the fin centerline tempera- 
tures from procedure (e). 

(g) Net heat flux to or from the passage walls by use of the integrated 
fluid temperatures and the wall centerline temperatures. 

(h) The values for metal thermal conductivity are obtained from the table 
at the passage wall centerline average temperature. 

The nature of the program makes it necessary to use a partly graphical 
solution for finding the relation between flow rate and fin geometry (weight) 
for panels of equal length. The general procedure used to determine flow rate 
as a function of fin geometry was the same for both high- and low-flux folded- 
flow systems. First, the intermediate hydrogen temperature between folded-flow 
and single-pass panels was determined by graphical techniques from the computer 
output. Then, the fin geometries that had acceptable maximum design wall tem- 
peratures and inlet pressures were determined, also by graphical techniques. 
Finally, plotting of effective fin thickness, t, as a function of flow rate for 
both panel systems permitted calculation of the cooling efficiencies. 

The intermediate hydrogen temperature between single-pass and folded-flow 
panels was calculated by the following procedure: From the computer output 

for various fin geometries in the folded-flow panel, maximum wall temperature 
was plotted as a function of panel length with parameters of hydrogen flow rate 
and outlet temperature (fig. 45). Cross plots of hydrogen flow rate and outlet 
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temperature, both versus maximum wall temperature, "® r ® f '? 

geometry as a parameter for a panel length of 2 ft (0.61 m)(fig. 46). Also 
from the computer output for single-pass flow, plots of hydrogen temperature as 
a function of length were prepared for various fin geometries with 
a parameter (fig. 47). The single-pass panel hydrogen inlet temperature wa 
then read as a function of flow rate at 2 ft (0.61 m) upstream from the length 
at which the design value of hydrogen outlet temperature occurred. Hydrogen 
outlet temperature from the folded panel and hydrogen met temperature to the 
single-pass panel were then plotted in fig. 48 for low-flux panels and fig. 49 
foAigh-flux panels. The intersection of these curves is the locus of inter- 
mediate hydrogen temperatures for the flow rate and fin geometry combinations 
indicated. The narrow range of hydrogen flow rates that satis Y a avail- 

conditions of hydrogen temperatures, heat flux, and fin geometry was then avail 

able. 

The procedure for determining configurations with acceptable maximum wall 
temperature or inlet pressure is described next. For low-flux, folded-flow 
panels, the maximum wall temperature was plotted vs flow rate with 1 " V 
as a parameter (fig. 50) to find those geometries that had maximum wall temp 
tures below I860°R (I033°K). This maximum wall temperature was based on the 
use of Concept I (fig. 51 ), where the heat exchanger and structural panel are 
combined into a single layer, and the maximum design value for cross section AT 
was 1 00° R ( 55. 5° K) . Reference to fig. 50 shows that folded-flow heat exchangers 
with 14 fins/in. (5.5 fins/cm) or less had acceptable maximum wall temperatures. 
The maximum wall temperature in the single-pass panel was less than I860 R 
( 1 033° K) for all fin designs in fig. 48. 

The maximum pressure drop for all low-flux panel systems was 2 P^' ( > 3. 8 
kN/m 2 ) with an outlet pressure of 300 psi (2070 kN/m 2 ). The average tempera- 
ture used for pressure drop calculation in the panel system with a folded-flow 
and a stngle-pass panel was obtained from the inlet and outlet temperatures to 
three flow lengths, each of 2 ft (0.61 m) . 


Average T c = (T CJ + 2T HJ + 2T HM + T CQ )/6 


( 20 ) 


For the single-pass panels in parallel, the average temperature was obtained 
from: 


Average T c = (T^ + T CQ )/2 


Inlet pressure for the 
series at high heat flux was 


system with a folded and a single-pass 
found from eq. 21 y as follows: 


panel 


i n 


P CI = p co + AP c = 


CO 


+ 11 + 


^) (R, ( 


T CI + 2 T HT + 2 T HM +T CO p A p 




0 . 5 

( 21 ) 


where the terms are defined in the nomenclature except a&P i 
pressure drop and density ratio for the entire 6-ft ( I . 83-mj 


s the product of 
hydrogen flow 
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length. About half of the overall a AP occurs in the folded-flow panel with a 
flow length of 4 ft (1.22 m) and plain fins, and half occurs in the single-pass 
with a flow length of 2 ft (0.61 m) and offset fins. 

Inlet pressure for the single-pass panels was calculated from eq. 22: 

0 . 5 

( 22 ) 


p ci p co + AP c 


CO 


I + 


AP 

■ 

AP 


7 > 


<r) < t ci + V 




(*) 
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APPENDIX B 


APPLICATION TO SPECIFIC CONFIGURATIONS 


Three basic concepts shown in fig. 51 were considered In depth over a 
range of environmental conditions in the tradeoff studies of ref. I, Ranges of 

heat flux and pressure used as panel design points are summarized along with 

the resulting hydrogen flow and panel thickness and weight. Heat exchanger 

weight is noted to indicate the relative importance of minimizing this fraction 

of the total . 

To permit panel weight minimization during concept evaluation and tradeoff 
studies presented in ref. I, figs. 52 through 60 were prepared for single-pass 
panels. Either (l) design maximum wall temperature or (2) the difference 
between design maximum wall temperature and hydrogen outlet temperature is 
plotted vs fin height for a range of fin thicknesses or fins per inch. Where 
pressure containment could not be achieved by fins of minimum gage and spacing, 
the plots of hydrogen inlet pressure and fin allowable pressure were used. The 
heat transfer aspects for obtaining these curves are discussed here. 

The boundary conditions established during the tradeoff studies and based 
on structural analysis include those listed in the Statement of Problem. A 
heat exchanger hydrogen outlet pressure of 300 psia (2070 kN/m 2 ) was used to 
allow a 50 psi (345 kN/m 2 ) outlet manifold pressure drop. A heat exchanger hot 
wall thickness of 0.010 in. (0.0254 cm) and thermal conductivity of 14.5 Btu/hr- 
F-ft (25.1 W/m-°K) were used with the hydrogen outlet temperatures, dimensions, 
concepts, and heat fluxes noted in figs. 52 through 60. The maximum tempera- 
tures of the adiabatic structural panel were assumed equal to the selected 
values of hydrogen outlet temperatures. 

The design maximum wall temperature was usually below the maximum value of 
2000°R (ill l°K), especially at the lower end of the heat flux and hydrogen out- 
let temperature range. If the maximum structural temperature of the hydrogen 
outlet had been allowed to increase and the design maximum wall temperature had 
been held constant, the hydrogen flow rate would have reduced from the values 
indicated, but the structure weight would have increased above the values noted 
i n ref. I . 

Rectangular offset fins and plain round tubes were the heat exchanger 
geometries studied for Concepts 2 and 3, fig. 51. The geometries for Concept I 
are plain rectangular fins. 

Inlet pressure for Concepts I and 2 was calculated by: 



where crAP/A may be read at W/w from figs. 12 through 15 for all concepts. 
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Inlet pressure for Concept 3 was calculated by: 


P CI = p co + AP c = 


l 

AP \ 

D 2 . I 

1 + — 1 

P C0 + \ 

APf / 

_ 



( R )(o.9t ci + ugi^jd) 


'W 


0 . 5 


(24) 


As indicated by eq. (24), the pressure drop for Concept 3 is always 
slightly higher than for Concept 2 because the hydrogen inlet temperature. to e 

hot heat exchanger is higher by 10 percent of the hydrogen AT, thus 
density. This AT represents a nominal heat leak to the cold ^tructura P 
tion heat exchanger. A value of 1.06 was used for the parameter (I + AP m /AP f ) 
in all calculations of P^. 

For Concept 3, the pressure drop in the cold structural protection heat 
exchanger was made negl i gi ble.( 1 ess than I percent of hot-heat exchange r^P)^ V 
use of a plain aluminum fin with a height of 0.050 in. (0.1 )• 

no indentive to design a minimum weight fin for the cold heat exchanger hecaus 
elimination of the fin altogether produces a we.ght reduct. ° n w ( j i be 

lb/ft 2 (0.34 kg/m 2 ). The hydrogen pressure in the cold heat exchange r wi 
near the design maximum of 1000 psia (6900 kN/m') because the lightest we.ght 
h ", exchangers for high-flux panels are generally ng th« ( h. gh«t 

permissible inlet hydrogen pressure and pressure r°P- ’ exchanqer 

selected for hydrogen pressure contai nment strength in the coldheat ^changer 
is 20(7. 9)R - 0.050(0.127) - plain - 0.005(0.013) - (7 - 3, L Q /4r h , 

constant y . This combination of fin spacing and fin thickness provides the 
metal cross section required for pressure containment. 


Design maximum wall temperature for all concepts was calculated by: 


T DMW T C0 + AT fin + 3 ^ T WH T W^ 


q/A 


2q/A 


(25) 


CO T^hAy/wA 3k 

where y^/w, may be calculated at T [0 by the method described in Appendix A 
or may be read at W/w from fig. I 2 through 15 and modified by: 

( T| hA-p/w-O at 2000°R (llll°K) from figs. 1 2 through 15 
(T) o hA T /w^) at T cq = — 2 “ C 

(m hA /w£) at 2000°R (llll°K) from fig. 21 
and c = (T) o hA T /wJi) at T CQ from fig. 21 


(26) 


, , ..a. cc , r _ is required only when all of the 
‘o’" T CO 


Use of fig. 21 to obtain Tl Q hAj/wA at T lu 

thermal conductance data are at one temperature as in fi 9 S- '2 through !5. If 
thermal conductance is calculated at the hydrogen temperature of ^terest th 
use of C is not required. Reasons for use of 2/3 (T yH - T y ) are 
ref. I. In ref. I, T was used to determine the fin material strength and fin 

allowable pressure of ° the various heat exchanger geometries considered. 
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The uniform heat fluxes shown on the curves for Concept 3 were used for the 
hot heat exchanger. As noted previously, 10 percent of this heat flux was 
assumed to pass through to the cold structural protection heat exchanger. Con- 
sequently, 10 percent of the hydrogen temperature change was assumed to occur 
in the cold structural protection heat exchanger. The remaining 90 percent of 
the hydrogen temperature change occurred in the hot heat exchanger. The heat 
transferred from the hot to the cold heat exchanger would cause an increase in 
cooling efficiency at less than infinite recovery temperatures. 

The hydrogen inlet pressure P is not usually indicated at a heat flux of 

less than 100 Btu/sec-ft 2 (1135 kW/m 2 ) because the hydrogen pressure drop is low 
enough to be practically negligible. For example, in Concepts 2 and 3 the 
hydrogen pressure drop at this heat flux was always less than 15 psi (103 kN/m 2 ) 
with an outlet pressure of 300 psia (2070 kN/m 2 ). The pressure drop is not 
reported for Concept I because the greatest pressure drop caused by any fin 
geometry selected for Concept I was 56 psi (386 kN/m 2 ). This pressure drop 
occurred for a length of 5 ft (1.53 m) with a fin geometry of 20(7. 9)R - 0.050 
(0.127) - 60(153) - 0.003(0.008) - (7 - 3). 

For the plain, round, tubular heat exchanger geometry, curves for two dif- 
ferent hydrogen flow rates are shown in figs. 57 and 58, even though only one 
panel length and heat flux were used for each figure. The 18-percent higher 
flow rate is based on the assumption that the heat flux may be increased by the 
ratio of total exposed area to panel projected area. This is commensurate with 
the heat exchanger geometry description given in the Symbols and Parameters 
table. A tube wall thickness of 0.010 in. (0.025 cm) was used to have the same 
hot wall thickness used for the corrugated-fin heat exchanger geometries. 

To compare the performance of bimetal and single metal fins, performance 
was calculated for rectangular offset fins with a thermal conductivity of 100 
Btu/hr-°R ft (173 W/m°K) and is reported in figs. 59 and 60 for the minimum 
weight geometries. The fin thickness was assumed to be made up of equal parts 
of copper and Hastelloy X. The pressure containment strength was based on the 
Hastelloy X portion of the fin only. While the higher fin conductivity pro- 
duced lower design maximum wall temperatures, the reduction was not enough to 
produce lighter weight heat exchangers than for a thermal conductivity of 10 
Btu/hr-°R-ft (17.3 W/m°K). The allowable fin pressure dictates the minimum 
weight design in fig. 59; the maximum hydrogen inlet pressure dictates minimum 
weight design in fig. 60. 

For performance evaluation where the recovery temperature is 5000°R 
(2780°K), the hot-gas heat transfer coefficients were calculated by eq. I. 

These heat transfer coefficients were used in the four-fluid heat exchanger 
analysis program, Appendix A, to calculate the data plotted in figs. 61 through 

64 for a length of 2 ft (0.61 m) . Hydrogen flow rate, cross-section AT, and 
pressure drop read from these figures at hydrogen outlet temperatures of 1400, 

1 600, 1760, and I900°R (778, 888, 976, and |055°K) were used to prepare figs. 

65 through 69, which are plots of design maximum wall temperature and hydrogen 
inlet pressure as functions of heat exchanger geometry. On the low-flux curves, 
where pressure drop is not significant, the hydrogen pressure containment 
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capability of each geometry is noted at one value of design maximum wall 
temperature. At the higher fluxes, where pressure drop is significant, the 
allowable fin pressure is plotted as a curve. 


of figs. 61 through 64 is that hydrogen flow 
in fin height and fin spacing only at fluxes 


A noteworthy characteristic 

rate is reduced with an increase - .... • 

above 10 Btu/sec-ft 2 (114 kW/m 2 ). Increased fin height and spacing cause lower 
thermal conductance, higher fin AT, and higher wall temperature along the entire 
thus reducing the heat flux into the panel. The hydrogen 


panel length, 
rates for the various 


geometries noted in figs. 61 through 68 for Tj. 


flow 
'of 5000° R 


( 2780°K) are always less than for the similar size and nominal heat flux panel 
with infinite recovery temperature. 


Off-Design Performance 


the fin 
In al 1 


Steady-state performance was calculated for minimum-weight-panel heat 
exchangers selected during the tradeoff studies at various design heat fluxes. 
The results of these calculations are presented in figs. 70 through 7 , 
geometry, design heat flux, and flow rate are noted on these figures, 
cases, the heat fluxes range between 10 and 1 00 percent of design heat lux 
since maximum thermal stresses occur at design heat flux. The mi mmum weight 
fins analyzed for Concept I are for a normal pressure differential of 6.95 psi 
(47.9 kN/m 2 ). The normal pressure has no effect on the other heat exchanger ^ 

geometry selections 


Constants for the heat exchanger geometries are hydrogen 
I600°R (888°K), flow length of 2 ft (0.61 m) , infinite 
and fin thermal conductivity of 10 Btu/hr-°R-ft (17.3 


outlet temperature of 
recovery temperature, 

W/m-° K) . 

The design maximum wall temperature and cross section AT curves for plain 
fins in figs. 72 and 73 do not increase uniformly as flux increases, but have a 
maximum value at a flux less than design point flux. The heat transfer per- 
formance curves from fig. 7-3 of ref. 2 used in estimating the plam-rectangular- 
fin performance for Concept I have a sharp variation at Reynolds numbers in the 
transition region between laminar and turbulent flow. Trans 1 ti on be gi ns at a 
Reynolds number corresponding to a heat flux of about 70 Btu/sec-ft 2 (796 kW/m ). 
There is no similar sharp transition and inflection i n the performance curve 
from fiq. 10-61 of ref. 2 used for rectangular offset fins in Concepts 2 and 3. 
The design maximum wall temperature and cross section AT. curves. are, therefore, 
smoothly increasing functions of heat flux. The uncertainties in per ormance a 
transition Reynolds numbers for plain geometries are a compelling reason for 
exercising caution when designing for their use. Flow instability is a potential 
problem at high ratios of outlet-to-i nlet temperature, and at laminar Reynolds 
numbers, especially in plain geometries. 
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APPENDIX C 


METHOD OF MANIFOLD ANALYSIS 


The flow rate parameter used for data presentation in connection with the 
manifold analysis combines heat flux and panel dimensions. The parameter is 
related to hydrogen flow rate, specific heat, and temperature change by eq. 27. 


( q/A) ( 1 / w) 


(W/")(C P KT C0 - T C1 ) 

w 


(27) 


Eq. 27 was obtained by rearranging and dividing eq. 2, Appendix A, by w. The 
parameter ( q/A) ( i,/w) provides a method of indicating flow rate without requiring 
a separate curve for each flow rate in each panel width. Inspection of eq. 27 
shows that a separate flow rate is associated with each curve of the type shown 
in fig. 38; but each flow rate applies to a broad range of combinations of panel 
heat flux and dimensions. 


To calculate pressure drops, the manifold resistances were separated as 
shown in fig. 76. The pressure drop In the separate resistances was calculated 
by application or extrapolation of loss coefficient data from ref. 6 for incom- 
pressible flow in mitered bends with expansion or contraction; and from directly 
applicable data in ref. 2 for friction in plane rectangular channels. The pres- 
sure drops and the separate resistances were added and the results are plotted 
in fig. 77 for either inlet or outlet manifolds in three wi dth-and-port-di ameter 
combinations and for two manifold fin heights. The pressure drop data in fig. 

77 is for any combination of heat flux and 1 ength-to-wi dth ratio for fixed mani- 
fold widths. If a manifold is to be designed for a fixed heat flux and length- 
to-wi dth ratio, then a cross plot of fig. 77 can be made as in fig. 78, which 
applies to a broad range of manifold widths at discrete values of the flow rate 
parameter, ( q/A) ( jJ/w) . 

The pressure drop for both manifolds combined was calculated from eq. 28 
and plotted in fig. 38, as well as in fig. 79, for a broad range of manifold 
widths and values (q/A)(£/w). 


AP 


M 


jAP MI , J&P H0 
P MI P M0 


(28) 


The important geometric and hydrogen pressure and temperature limitations 
are noted in fig. 70. Calculations were made for both inlet and outlet mani- 
folds but the resultant values for aAP^j and ctAP^q never differed by more than 

15 percent. Since a tolerance of ±15 percent on the calculated loss coefficients 
for the various resistances is realistic, the larger value was assumed to apply 
for both inlet and outlet manifolds. 


Table 7 indicates the fraction of manifold pressure drop contributed by the 
various resistances in fig. 76 for three manifold widths. The bend-wi th-area- 
change and horizontal friction caused from 40 to 90 percent of the pressure drop 
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in both the inlet and outlet manifolds. The last three resistances in the 
outlet manifold caused from 81 to 94 percent of the pressure drop in both mani- 
folds combined, since the inlet manifold causes only 3 percent of the combined 
pressure drop. Also, from table 7, the outlet port pressure drop was calculated 
to be 40 percent of the overall manifold pressure drop. The remaining 60 per- 
cent of the overall manifold pressure drop is inversely proportional to the 
square of the manifold fin height. The manifold fin height controls the free 
flow area, A f , which in turn controls pressure drop according to eq. 29. Eq. 29 

applies to subsonic Incompressible flow as it exists in panel manifolds with 
reasonable pressure drop. 


= (Zk)W 

" A f 2 v 


(29) 


To calculate manifold fin height it was first necessary to determine mani- 
fold pressure drop as a function of the flow parameter (q/A)(A/w) for the 
specific manifold widths to be used in the tradeoff studies of ref. I. These 
manifold pressure drops are shown in fig. 36 for a hydrogen outlet pressure of 
250 psi (1560 kN/m 2 ). The manifold fin height was then calculated by using 
eq. 30. 


new 


i n M 


basi c h 


I 0.6 AP c 

f i n M y AP^ - 0. 4AP 


-r 


(30) 


where basic h r M = 0.25 in. (0.635 cm) 
tin M 

AP c = manifold pressure drop from fig. 36 

AP = pressure drop allotted to outlet manifold of 45 psi 
M (310 kN/m 2 ) 

While the results discussed above show manifold pressure drop for uniform 
flow, the flow distribution resulting from existing pressure drop distributions 
can be calculated by similar methods. The analysis is iterative as used here. 

In the analysis, uniform flow was assumed so that each flow path in the 
manifold from the longest at the edge to the shortest at the center received 
equal flow. The overall pressure drop of each local flow path, AP i q C a 1 J ^^ rou 9^ 

both manifolds and the core was calculated by eq. 29 and used to obtain a first 
trial flow-distribution profile, related to the shortest path by use of eq. 31. 

AP, ,\ 0 - 5 /AP, , 

local \ / local 

AP rn / I AP • 

SP / \ mi n 

To obtain eq. 31, the terms are eliminated from eq. 29 that are constant for any 

specific manifold and core geometry with any conditions of hydrogen flow rate, 
pressures, and temperatures. By use of the nonuniform flow profile obtained 


W 


1 oca 1 
W 

max 


W 


local 


W 


SP 



48 


from eq. 31 , the pressure drop through all passages from the longest at the 
edge, AP Lp to the shortest at the center, AP $p Is calculated. A second non- 

uniform flow profile was then obtained by use of eq. 31 with the new pressure 
drop profile. This iterative procedure is continued until the same pressure 
drop is calculated for all flow paths. This process is based on the principle 
that all flow paths in parallel between two common plenums or pipes have the 
same pressure drop. The flow distribution profiles resulting from the above 
iterative calculation procedure are shown in fig. 39 for the three manifold 
widths. The set of typical operating conditions and manifold geometry shown 
in table 5 was selected to allow preparation of fig. 39. 



APPENOIX D 


MANIFOLD TESTS 


This appendix describes the test apparatus and procedures used to evaluate 
manifold flow distribution and pressure drop. Data reduction methods, data 
presentation, and detailed test results are discussed. 


Test Apparatus 


The test unit,, shown disassembled in fig. 80 and installed with instrumen- 
tation.^ fig. 81, consists of inlet and outlet plenums, inserts for two mani- 
fold fin heights, and a 65-tube core. The core is split to allow Insertion of 
an orifice plate flow resistance. The Installation of the orifice plate is 
optional, and thus permits operation at two levels of core flow resistance. 


Overall inlet and outlet pressures are obtained by instrumenting each 
plenum with five static pressure taps. The core contains three rows of 16 
static pressure taps on individual tubes. The larger (Inlet) core section con- 
tains two rows of pressure taps from which are obtained core Inlet pressure and 
pressure drop data. The smaller core section contains the core outlet row of 
pressure taps. Fig. 82 schematically depicts the Instrumented test unit with 
the core geometry noted. 


The readout of the static pressures at the inlet and outlet plenums, and 
the core inlet and outlet (Stations 2 and 4) is accomplished on a 48-port 
Scanni valve. This valve may be indexed to each port location; the indexing 
feature provides pressure correspondence between the center tap and the selected 
port. Station 3 pressures are fed to a common manifold and individually valved 
off by means of toggle valves. Core pressure differences were measured by 
selecting the appropriate Scannivalve and Station 3 locations. A water manometer 
board was used for core pressure drop readouts between Stations 2 and 3, and 
pressure gages and a mercury manometer board are used for all static pressure 
readouts . 


The supporting apparatus consists of flow control, filtering, and measuring 
sections upstream of the test unit. Flow measurement is accomplished upstream 
of the test unit by means of a sharp-edged orifice. The isothermal airflow tem- 
perature is recorded at this measuring section. 

Photographs of typical manifold test specimens are presented in figs. 83 
and 84. Specimen geometries, as measured from the test specimens, are given in 
table 6. The flow distribution insert Is described later in this section. The 
sequences of parameters that were compared are summarized in the following table. 
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Comparison parameter sequence 
Port diameter 

Port-radi us-to-di ameter ratio 
Man i fol d wi dth 

Fin height and geometric similarity 
Rectangular and tapered manifold comparison 
Flow distribution insert 


Specimen numbers 
5, 2, and 4 
7, 6, and 2 
3, 2, and I 
12, 2, and I 
II, 10, 2, and I 
I with insert 


The manifold flow passages are formed with mitered, plain rectangular fins, 
sandwiched between two manifold plates. The fins would be brazed between thin 
sheets about 0.010 In. (0.025 cm) thick to provide adequate coolant pressure 
containment with light weight in a flight hardware application. Therefore, only 
Specimens I and I I were tested without fins for reference. The fins were 
attached to one manifold plate with a rubber-base cement and removed as neces- 
sary for tests by soaking in toluene. Balsa wood inserts at the sides and 
bottom of the manifold were used between plates to prevent bypass leakage when 
fin configurations were bolted into the test unit. 


Data Reduction Method 

Flow distribution across the core was determined from the measured static 
pressure drop distribution between Stations 2 and 3. The basic equation relat 
ing measured pressure drop and local flow rate in any tube is: 



where the terms are defined in the table Symbols and Parameters and 
A = 0.558 in 2 (3.6 cm 2 ) 


l = 5 in. ( I 2. 7 cm) 


r = 0.1045/4 = 0.0261 In. (0.0662 cm) 
h 

p = P 2 + P 3 /R(T 2 + t 3 ) 

H a 2 - 3 


The core friction factor, f, 
uniform ai rflow. The flow rate i 


is calculated from core test results with a 
n each of 16 tubes with pressure taps is cal- 


w 


1 ocal 


2 9p AP 2-3 A 2 

a 2~5 


0.5' 


2p t-L - ±\ ♦ s 

a *-3 \ps 111) r h 


The flow distribution is plotted as 


^local ^ ^local 


W 


W A 


avg ' tota 1 
where 65 is the number of core tubes. 


(33) 


(34) 


Eqs. 33 and 34 form the basis of a computer program that was used to reduce 
the flow distribution test data. The required input includes core geometric 
parameters^ test static pressures and temperatures; total flow rate; and baro- 
metric pressure. Tables of core Fanning friction factor obtained from test unit 
calibration and tables of fluid properties are also included. The program 
prints out the core weight/flow distribution; as well as the distribution of 
associated parameters such as Mach number; Reynolds number; and total pressure 
and pressure drop. The arithmetic average core inlet and outlet total pressures 
are calculated from the core inlet and outlet static pressures and flow distri- 
bution. Using these average pressures and the test unit inlet and outlet total 
pressures calculated from test static pressures and flow, the program calculates 
the various total pressure drop parameters and entrance and exit loss coeffi- 
cients. In addition; the program plots the flow and pressure distributions. 


By choosing the plain tube core length between Stations 2 and 3 as the flow 
measuring section; the test friction factor calibration can be referenced to 
standard plain tube friction data (ref. 2). This referencing makes possible 
more accurate calibration data interpretation and extrapolation. Performance 
characteristics such as inlet and outlet manifold loss coefficients are deter- 
mined by using Stations I, 2; 4; and 5. 


Test Resul ts 

Test unit cal i brat ? on . --The tube Fanning friction factor is shown in fig. 
85. This characteristic was obtained by calibrating the core with nominally 
uniform flow. The corresponding flow distribution profiles are shown in fig. 86 
for various flow rates; with and without orifices in the tubes. The slightly 
nonuniform flow profiles (less than 5 percent) appear to be due to measurement 
scatter rather than to actual flow nonuniformity; since tests with a reversed 
flow direction also caused a reversal in flow profile tendencies. The flow 
profiles shown are for full core width. The left half of the core; where more 
uniform flow is indicated and where II of each group of 16 static pressure taps 
are concentrated; was used to calculate the flow distribution reported in figs. 

41 and 42, as well as elsewhere in this report. This use of test results from 
the left half of the core is also justified because the test unit and all speci- 
mens are symmetrical about the centerline. 
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The test points in fig. 85 include data from both the smooth tube section 
of 5-in. length and 9-in. overall length without the 0. 063-i n. -di ameter orifice 
plate. Data from ref. 2 for 1 ength-to-di ameter ratios, l/k r^, of 100 and 

infinity are plotted to assist in correlating the data. The actual 1 ength-to- 
di ameter ratio between Stations 2 and 3 is 48, but the effective 4/4^ is 

between 1 00 and infinity, since the tube entrance is 19 diameters upstream of 
Station 2. Agreement is good throughout the turbulent Reynolds number range. 

The test data near the transition Reynolds number range show a slight upward 
curvature not present in the reference curves. 

Core flow prof i les . --As mentioned above, only the test data from the left- 
hand half of the core were reduced (with the exception of the test unit calibra- 
tion data). Fig. 87 is representative of the static pressure data obtained for 
half the core and used with the data reduction method described above to obtain 
the flow profile for Specimen I (see fig. 89). The 0.0 core flow width coordi- 
nate is at the edge of the core and the 1.0 coordinate is at the centerline of 
the core. 

In general, the tests for which results are presented in fig. 88 through 
93 have identical manifold specimens at the core inlet and outlet. For addi- 
tional information, some flow distribution curves in fig. 92 are presented for 
Specimen I at the inlet only and at the outlet only; and individual Specimen I 
samples have been designated as Specimen IA and IB to indicate the effects of 
manufacturing tolerance. These figures include results from tests with the flow 
distribution insert. Design and test of the insert are discussed later in this 
section. As expected, tests on Specimens I and II with the manifold fin removed 
do not exhibit the flow peak at the core centerline present in tests with fins, 
as shown in fig. 90. Test results with Specimen I at the inlet and Specimen 2 
at the outlet are similar to those obtained with Specimen I at both inlet and 
outlet. 

Unless otherwise noted, all reported tests are with the core orifice 
resistance; all manifold loss coefficients are based on the core area, test 

manifold loss coefficients, K and K Q , are based on uniform flow in the core 

area and test unit inlet or outlet density (Stations I and 5, respectively). 

The flow distribution is expressed as a ratio of local flow in the core to 

average flow in the core. The flow profiles in figs. 89 through 93 are at a 

nominal airflow rate of 4 lb/min (0.0303 kg/sec). Fig. 88 shows a set of repre- 
sentative flow profiles for Specimen 2 at various flow rates. The flow profiles 
are almost identical because the flow in all tubes is in the turbulent Reynolds 
number regime. The minimum Reynolds number in any tube during any test was 
about 6000, just below the lower edge of the turbulent regime. 

Comparison of core and manifold flow prof i les . --General data correl ation 
is difficult for the wide range of geometric manifold parameters investigated, 
especially because of the hard-to-def i ne relationship between core and manifold 
flow profiles. The manifold configuration (table 6) provides opportunity for 
crossflow along the channel formed at the right-angle joint between the core 
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and the manifold. The flow at this joint adjusts itself to accommodate the 
static pressure profile along the joint. Because of the existing inequality 
in resistance of manifold paths, this pressure profile decreases from the 
location of the shortest manifold path to that of the longest path in the 
inlet manifold and from the location of the longest path to that of the 
shortest path in the outlet manifold. The crossflow is away from the short- 
est path in the inlet manifold and toward the shortest path in the outlet 
manifold. This crossflow pattern causes both inlet and outlet manifold flow 
profiles to be less uniform than those in the core. Calculations show, for 
example, that the manifold centerline flow is twice the core centerline flow 
for inlet manifold Specimens I and II. 

The core static pressure profiles are flat near the edge of the core for 
all specimens, which indicates that the crossflow there is small. Thus, the 
longest-path test loss coefficients reported in fig. 94 were evaluated assum- 
ing no crossflow at the core edge. 

A knowledge of the centerline manifold flow for inlet Specimen I permits 
analysis of the shortest-path pressure losses. In the longest path, the high- 
pressure-loss components were found to be the horizontal fin friction and the 
hori zontal -to-verti cal fin miter joint. Neither of these high-loss components 
is present in the shortest path. Calculations showed that vertical fin fric- 
tion and core entrance losses were not more than 15 percent of the measured 
total losses. The high pressure difference between the plenum and the core 
tube in line with the port centerline can be accounted for by postulating a 
nonuniform flow distribution in the inlet port; that is, flow through the inlet 
port adjacent to the shortest path is 3.5 times greater than the average port 
flow. 


General data correl ati ons . --A summary of test results for the various 
specimen configurations appears in table 8 and was used to prepare the data 
presentations. Figs. 41 and 42 are a summary of the experimental investiga- 
tion. The curves are drawn to indicate trends using easily identifiable param- 
eters rather than exact functional relationship. The general trend shown is 
that flow becomes more uniform as manifold pressure drop (resistance) is 

reduced. The maxi mum- to-mi ni mum flow ratio (W /W . ) and core-to-ove ral 1 

max min 

total pressure drop ratio ( AP2-4/AP1 -5) are, respectively, the characteristic 
flow and flow resistance uniformity parameters. Since the shortest-to-longest 
manifold path flow resistance ratio is small, the manifold pressure drop is 
rep resentat I ve of the nonuniformity existing in the various overall test unit 
flow paths. Hence, the related parameter, AP 2 -4/AP 1 - 5 * is a measure of the 
uniformity of flow resistance in the various paths. 

Calculated and test loss coef f i ci ents * --I t is helpful to know loss coeffi- 
cients for various parts of the manifold because flow distribution can be 
calculated if loss coefficients are known. Loss coefficients were calculated 
for comparison with test loss coefficients and to provide generally useful data 
for design of similar manifolds. In summary, calculated coefficients for the 
manifold longest path were from 0.87 to 1.84 times test coefficients (after 
friction was subtracted) for outlet manifolds and from 1.79 to 2.52 times test 
coefficients for inlet manifolds (fig. 94). 
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An indication of the relationships between calculated and test loss 
coefficients for the longest manifold path is shown in fig. 94. Both test 
and calculated coefficients are evaluated for the longest manifold path and 
are directly comparable. Test-effective loss coefficients (discussed in later 
paragraphs and presented in figs. 95 and 96) are related to the particular flow 
nonuniformity and, therefore, will have values that. i ncrease with increasing 
flow nonuniformity. For convenience in this analysis, the loss coefficients 
are based on the horizontal fin flow area because the majority of the longest 
path losses are related to this fin. The frictional loss characteristics 
available in ref. 2 for the horizontal fin are directly applicable and have 
been subtracted from the overall coefficient. The test coefficient was 
ated by using the measured local path pressure drop and airflow, assuming core 
and manifold longest path flow to be equal. 

Since the resistance of the shortest manifold path is small. in relation 
to that of the longest path, the resistance of the longest path is charac er- 
istic of the effective resistance of the manifold.. A theoret i cal cal cu at. on 
of this longest-path resistance and its relationship to the experimen a y 
determined effective manifold resistance can be used in the prediction of 
manifold pressure drop and, hence, of flow uniformity. 

The calculated longest-path loss coefficients and pressure drops are 
presented in table 9. The pressure drops are cal culated .on the 1 bas i s o a 
uniform 4 lb/min (0.0303 kg/sec) airflow through the manifold at a density 
of 0.0765 lb/ft 3 (1.225 kg/m 3 ). Also tabulated are the f 0 0303 

effective loss coefficients with a nonuniform test flow of 4 lb/min ( 0.0303 
kg/sec). Refs. 6 and 2 were used to estimate the longest-path pressure drop 
by dividing the path into the pressure drop components shown in table . 
indicated in this table, the effect of some of the longest-path components was 
aDDroximated by the sum of two or more pressure drop components (loss coeffi 
cien?s) w^en '"direct calculation of a single loss coefficient was unavailable. 

The relationship between calculated and test loss coefficients is P re " 
sented graphically in fig. 96. As previously mentioned, f ig . 9!> P ^®"* 5 
experimental relationship between the test loss coefficient and the flow ratio. 
Even though the curves drawn are not considered definitive, the trends a 
i ndi cated and, with the data points, can be used to predict the manifold losses 
and the flow nonuniformity of similar manifolds. 

The theoretical manifold flow resistance correlation presented in fig. 
shows the calculated longest-path resistance to be higher than the test- 
effective resistance. This tendency is qualitatively correct since the . 
lonqest-path resistance is always greater than the effective average resis- 
tance. Also, the calculated longest-path resistance diverges farther from 
the test-effective resistance as the resistances increase. ecause e 
shortest-path resistance is about the same for all specimens, an increasing 
manifold Resistance means an increasing longest-path to effect' ve-ave rage 
resistance ratio, and thus a diverging resistance curve. This trend 
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qualitative agreement with fig. 95, where an increased manifold resistance is 
shown to be accompanied by an increasing flow nonuniformity. 

Comparison of tapered and rectangular mani fol ds . — Fi gs. 90 and 9| show a 
comparison between f i nned-tapered and rectangular manifold flow profiles. A 
comparison of the associated manifold flow resistances is shown in fig. 95. 
Tapered and rectangular manifolds, Specimens II and I, respectively, are simi- 
lar in flow profile and flow resistance. This similarity is expected since 
the taper in Specimen II is small, and it approaches Specimen I in flow geom- 
etry. Tapered and rectangular manifolds, Specimens 10 and 2, respectively, 
however, exhibit considerable differences in flow profile and flow resistance. 
These differences are primarily a result of a narrower vertical fin flow path 
(between edges of the horizontal fin at the port) in Specimen 10 than in Speci- 
men 2 (fig. 83). This increased core centerline (shortest path) flow resis- 
tance in Specimen 10 causes a greater flow through longer paths, Increasing 
both pressure drop and flow uniformity. The same shortest-path resistance 
was assumed for all specimens while calculating the manifold pressure drop 
correlations in figs. 95 and 96. Hence, differences in shortest-path resis- 
tance, as emphasized by Specimen 10, cause data correlation scatter in these 
f i gures. 

Uniform flow manifold des i gn . --Uni form flow In the manifold can be 
obtained by equalizing the resistance through all manifold flow paths with an 
orifice plate inserted adjacent to the edge of the vertical manifold fin near 
the core. The holes in the orifice plate must be nonuni formly spaced or have 
various diameters, or both. Equal resistance in all paths ensures uniform 
manifold flow and pressure drop, independent of the core and other uniform 
seri es res i stances. 

A uniform-flow manifold design was based on test data for Specimen j, a 
rectangular manifold with one port. Nonuniformity in flow path resistances 
I s caused by 

(a) A variation in flow length (especially in the horizontal fin) 

( b) The absence of a right-angle mitered fin turn in the shortest 
manifold flow paths at the middle of the manifold 

The insert design shown in fig. 84 was attached at the exposed edge of 
the vertical fin. For a uniformly distributed 4 lb/min (0.0303 kg/sec) flow, 
the required nonuniformity in orifice plate pressure drop, as calculated for 
Specimen I, is shown in fig. 97. At the centerline, the flow paths have no 
right-angle bend, and this results in a step in the required pressure drop 
profile. The required pressure drops for both inlet and outlet manifolds are 
identical with the exception of those required near the core centerline. The 
difference at this location Is due to the higher calculated loss encountered 
in flowing through the expanding bend at the inlet port, as compared with the 
contracting bend at the outlet port. The insert tested was the same for both 
inlet and outlet manifolds. 
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Fig. 98 shows the calculated pressure loss characteristics of the orifice 
plate. The required orifice area is obtained by using figs. 97 and 98 and is 
shown in fig. 99. Using 0.030-in. (0.076-cm) diameter holes with a minimum 
clearance between holes of 0.020 in. (0.051 cm), the maximum avai 1 abl e area is 
0.014 in 2 per inch of orifice strip (0.029 cm 2 per cm of orifice strip). To 
obtain the larger required areas, the ends of the orifice strip are cut to a 
taper. A linear taper approximating the actual area variation is used. The 
hole density was maintained at greater than 10 holes per inch by using 0.020-in 
(0.051 -cm) diameter holes in the more restrictive locations near the port. 


Uniform Flow Manifold Tests 

The test results for Specimen I with the flow distribution insert are 
shown in figs. 92 and 93. The flow was not uniform, although the maximum- 
to-minimum flow ratio was reduced from 4.5 without inserts to 2.8 with 
inserts. The significance of the nomenclature used to identify the various 
curves is as follows: 

(a) Fig. 92 

(1) Inlet with outlet: Flow profile obtained with the basic mani- 

fold configuration of Specimen I, without flow correcting 

i nserts. 

(2) No. IA inlet: Flow profile obtained with the basic manifold 

configuration of Specimen I, used at the test unit inlet only. 
The test unit outlet is formed by a plenum without fins or 
insert and has a uniform resistance. 

(3) No. IB inlet: Same as for No. IA inlet, using the outlet mani- 

fold from (l), above, in the inlet for comparison purposes. 

(4) No. IA outlet: Flow profile obtained with the basic manifold 

configuration of Specimen I, from (2), above, used at the out- 
let only. The inlet is formed by a plenum without fins or 
insert and has a uniform resistance. 

(5) No. IA inlet with insert: The configuration described in (2), 

above, with insert added. 

(b) Fig. 93 

(1) Bottom curve: For basic manifold configuration of Specimen I, 

used in both inlet and outlet, with flow correcting inserts 
added to each, shown for comparison. Differences in the curve 
are due to manufacturing variances. 

(2) Top curve: For basic manifold configuration of Specimen I, 

with designated manifold ends run separately, the other end 
being at a uniform resistance. 
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Detailed modification of hole spacing can produce uniform flow. The 
insert was modified to improve the flow profile at the core center by enlarging 
the flow area in the center of the insert. The diameter of the seven center 
holes (approximately 0.5-in. (1.27-cm) span) was increased from 0.020 in. 

(0.051 cm) to 0.025 in. (0.0635 cm). Fig. 92 shows the flow distribution 
obtained. The max i mum-to-mi n i mum flow ratio was reduced from 2.8 with the 
basic insert design to 2.6 with the redrilled insert. The local change in 
profile, which was the purpose of the modification, was much more substantial. 
It does not appreciably affect the overall flow ratio because the maximum and 
minimum flows occur elsewhere. Further iteration in hole dimensions could be 
used to produce completely uniform flow in any manifold. 


59 



REFERENCES 


1. FI ieder , W. G.; Richard, C. E.; Buchmann, 0. A.; and Walters, F. M. : An 

Analytical Study of Hydrogen-Cooled Panels for Application to Hypersonic 
Aircraft. NASA CR-1650, 1970. 

2. Kays, W. M. ; and London, A. L. : Compact Heat Exchangers. Second Ed., 

McGraw-Hill Book Co., Inc. (New York), 1964. 

3. Anonymous: Properties of Parahydrogen. Report 9050-6S, Aerojet-General 

Corp. , PRA-SA-DSA-30, September 1963. 

4. Sparrow, E. M. ; Lloyd, J. R.; and Hixon, C. W. : Experiments on Turbulent 

Heat Transfer in an Asymmetrically Heated Rectangular Duct. ASME Paper 
No. 65-HT-8 . 

5. Elam, B. F.: Heat Transfer in Honeycomb-Core Sandwich Panels. ASME 

Paper No. 65-HT-I3. 

6. Aero-Space Applied Thermodynamics Manual. Society of Automotive Engineers 
Committee A-9, January 1962. 

7. Hendricks, R. C. ; Graham, R. W.; Hsu, Y. Y. ; and Freidman, R. : 

Experimental Heat Transfer Results for Cryogenic Hydrogen Flowing in 
Tubes at Subcritical and Supercritical Pressures to 800 psia. NASA TN 
D-3095. 

8. Hendricks, R. C.; Graham, R. W.; Hsu, Y. Y. ; and Medeiros, A. A. : 
Correlation of Hydrogen Heat Transfer in Boiling and Supercritical 
Pressure States. ARS Journal, February 1962, pp. 244-252. 


61 



TABLE I 


HEAT EXCHANGER GEOMETRY LIMITS 


Component 

Govern i ng 
cond i t ion 

Affected 

parameter 

L imi tat ion 

Heat transfer 
f i n 

Eros ion of 
parent meta 1 
by braze 
filler al loy 

Fin thickness 

Minimum thickness: 
superalloy: .003 in. (.008 cm) 

aluminum: .004 in. (.010 cm) 


Col lapse of 
fins due to 
braze 
f ixture 
loads 

Fin height, 
density, and 
th i ckness 

Minimum collapsing load 
(34.5 kN/m 2 ) at 2I00°F 
for superal loys 

of 5 ps i 
( I420°K) 


Formi ng 
tools and 
mater i a 1 
properties 

Fin density 
and 

th i ckness 

Max. fin 

Fins/in. thickness 

(Fins/cm) in. (cm) 

40 (15.8) max. .003 ( .008) 

30 (11.8) .004 (.010) 

20 ( 7.9) .006 (.016) 

10 ( 3.94) .010 (.025) 


Fabr i cat ion, 
pressure drop 

Fin he i ght 

.025 in. (.063 cm) min. 



Weight, 
fin AT 


.100 in. (.254 cm) max. 


Tubes 

Fabr i cat i on 
(assembly 
and brazing) 

Outs i de 
d i ameter 

.050 in. (.127 cm) min. 
.100 in. (.254 cm) max. 


Heat exchanger 
surface sheet 
(tube wall) 

Hand 1 i ng, 
particle 
damage, 
fabr icat ion 

Th i ckness 

Mini mum th i ckness : 
superalloy: .010 in. (.025 cm) 

aluminum: .016 in. (.040 cm) 

Coolant 
man i fo 1 d 
pi pes 

Hand 1 i ng 

Wall th i ckness 

Mi ni mum thi ckness 
. 030 in. ( . 0762 cm) 



63 







TABLE 2 


COOLING EFFICIENCY FOR FLOW ROUTING CONCEPTS 


Heat flux, 
Btu/sec-f t 2 
(kW/m 2 ) 

Recovery 

temperature, 

°R 

(°K) 

Cooling efficiency 

Single 

pass 

Folded- in- 
depth 

Fo 1 ded- i n- 
wi dth 

10 

3000 

0.47 

0.69 

0.64 

(114) 

(1670) 





5000 

.78 



.84 


(2780) 





7000 

.83 

.89 

.94 


(3890) 




250 

5000 

.81 

— 

.84 

(2840) 

(2780) 





7000 

.89 

.86 

.89 


(3890) 




500 

5000 

.84 

— 

.86 

(5680) 

(2780) 





7000 

.86 

.86 

.88 


(3890) 





64 



















TABLE 4 


MANIFOLD GEOMETRY FOR 
TRADEOFF STUDIES 


Flux, 

Btu/sec-f t 2 
(kW/m 2 ) 

Length, 

ft 

(m) 

(q/A)U/w) 
Btu/sec- f t 2 
(kW/m 2 ) 

No • of 
mani fold 
segments 
(ports) 

Manifold 
segment 
wi dth, i n. 
(cm) 

A V 

psi 

(kN/m 2 ) 


Port 

d i ameter , 
i n . 
(cm) 

250 

(2840) 

2 

(.61) 

250 

(2840) 

1 

24 

(61.0) 

20 

(138) 

. 142 
(.36) 

1.75 

(4.45) 

250 

(2840) 

3 

(.914) 

375 

(4260) 

1 

24 

(61.0) 

42 

(290) 

.236 

(.60) 

1 .75 
(4.45) 

500 

(5680) 

2 

(.61) 

1000 

( 1 1 350) 

2 

12 

(30.5) 

32 

(221) 

. 193 
(.49) 

1.25 
(3. 18) 

500 

(5680) 

3 

(.914) 

1500 
(17 030) 

2 

12 

(30.5) 

72 

(496) 

.407 

(1.03) 

1.25 

(3.18) 

500 

(5680) 

3 

(.914) 

2250 

(25 540) 

3 

8 

(20.3) 

55 

(380) 

.300 

(-76) 

1 .00 
(2.54) 
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TABLE 5 


MANIFOLD GEOMETRY AND OPERATING CONDITIONS FOR 
CALCULATED FLOW DISTRIBUTION 


Panel length, ft (m) 2 (0.61 ) 

Heat flux, Btu/sec-ft 2 (kW/m 2 ) 250 (2840) 

T CI' ° R 100 (55.5) 

T C0' ° R (° k ) 1600 (888) 

W/w, 1 b/sec-f t (kg/sec-m) 0.0935 (0. 139) 

Core pressure drop, ps i (kN/m 2 ) 40 (276) 

P CI , psia (kN/m 2 ) 500 (3450) 


Manifold geometry: 

Horizontal fin 

8 (3. 1 5 )R - 0.164 (0.417) - plain - 0.004 (O.OIO) 
Vertical fin 

10 (3.94)R - 0.164 (0.417) - plain - 0.004 (O.OIO) 



MANIFOLD SPECIMEN GEOMETRY 










TABLE 7 


MANIFOLD LOSS COMPONENTS 


Percent of overall manifold pressure drop 


Manifold width, ft fm) 

r- -* \ I 7 \ f / 


Loss components 

■HBUEH 

1 (.305) 

2 (.61) 

In 1 et man I f o 1 d : 




Inlet port 

14.0 

10. 1 

5.2 

Horizontal friction 

20. 1 

31.3 

48.4 

90-deg turn and expansion 

50.0 

50.9 

43.4 

Vertical friction 

2.2 

.9 

.3 

Three 90-deg turns 

9.6 

4.7 

1.9 

Contraction at core inlet 

4. 1 

2. 1 

.8 


100.0 

100.0 

100.0 

Outlet manifold: 




Expansion at core outlet 

4.4 

2.2 

.9 

Three 90-deg turns 

9.6 

4.7 

1.9 

Vertical friction 

2.2 

.9 

.3 

90-deg turn and contraction 

18.3 

16.9 

1 1.5 

Horizontal friction 

20. 1 

31.3 

48.4 

Exit port 

45.4 

44.0 

37.0 


100.0 

100.0 

100.0 


( q/A) ( l/vi) = 250 Btu/ f t 2 -sec (2840 kW/m z ) 

Fin geometry = 10 (3.94)R - .250 (.635)- plain - .004 (.010) 
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TABLE 9 

CALCULATED LONGEST-PATH AND TEST-EFFECTIVE MANIFOLD 
RESISTANCES AND LOSS COEFFICIENTS 


Specimen 

1 

2 

3 

4 

5 

6 

7 

10 

1 1 

ID 

Inlet manifold components 










■SiS 

A. Port 

K (port)* 

1.2 -e 





.75 

.39 






-p— 1 

1.2 

1.2 

1 .2 

AP (psl) 

9.3 

2.3 

t.o 

0 . 1 

6.5 

1.4 

0.8 

2.3 

9.5 

0.6 

(kN/m 2 ) 

64.2 

15.9 

6.9 

.69 

44.8 

9.7 

5.5 

15.9 

65.5 

4. 1 

8. Horizontal fin friction 











AP ( ps i ) 
(kN/m 2 

103.3 

713. 

12.9 

89. 

4. 1 
20.2 

4. 1 
28.2 

19.3 

133. 

6.4 

44.2 

4. 1 
28.2 

12.0 

82.8 

1 15. B 
799. 

4.6 

31.7 

C. Fin mi ter joint 










K (horizontal fin)* 

2.40 

2.22 

2.05 

1.93 

2.30 

2.07 

1.93 

2.04 

2.04 

2.22 

AP ( ps i ) 

95.3 

22.2 

9.9 

6.3 

34.4 

10.2 

6. | 

19.0 

90. 1 

9.0 

( kN/fn ) 

658. 

153. 

60.4 

43.4 

238. 

70.4 

42. } 

131 . 

621. 

62. 1 

D. Vertical fin friction 










AP (psi) 

0.4 

0.4 

0.4 

0.6 

0.4 

0.5 

0.6 

0.4 

0.4 

0.0 

(kN/m Z ) 

2-76 

2.76 

2.76 

4. 14 

2.76 

3.45 

4. 14 

2.76 

2.76 

0 

E Core entrance 











K (vert ical fin)* 

1.6 < 

3 









AP (psi) 

0.9 







— t 

:> 1 . 6 










0.3 

(kN/nr ) 

6.2K3 







<6.*, 












2.07 

Overa 1 i 











Calculated longest path 
AP (psi) 

209.2 

38.7 

16.3 

12.0 

61.5 

19.4 

12.5 

34.6 

2 16.7 

14.5 

(kN/m 2 ) 

1442. 

267. 

1 12. 

B2.B 

424. 

134. 

86.2 

238. 

1495. 

100- 

Kj* Calculated longest path 

500. 

92.4 

38.9 

28.8 

147.2 

46.5 

30.0 

82-6 

5 16. 

34.6 

Test effect i ve 

215. 

46.7 

25.4 

18.6 

91.4 

35.3 

2 1.9 

75.9 

278. 

18.2 

Outlet manifold components 











A . Core ex i t 











K (vertical fin)* 

1.6 < 

I 









AP (psi) 

0.9 







D 

>■1.6 











0.9 

0.3 

(kN/nr) 

6.2 1 







6.21 

2.07 









1> 

B. Vertical fin friction 











AP (psi) 

0,4 

0.4 

0.4 

0.6 

0.4 

0.5 

0.6 

0.4 

0.4 

0.0 

(kN/m 2 ) 

2.76 

2.76 

2.76 

4. 14 

2.76 

3.45 

4. U 

2.76 

2.76 

0 

C . Fin miter joint 











K (horizontal fin)* 

2.08 

2.06 

2.01 

1.97 

2.07 

2.02 

1.97 

1.88 

1.72 

2.06 

AP (psi) 

82.8 

20.6 

9.7 

6.4 

31.0 

10. i 

6.0 

17.7 

76.0 

8.3 

( kN/m 2 ) 

571. 

142 

66.9 

44. | 

214. 

69.6 

41.4 

122. 

524. 

57.2 

D. Horizontal fin friction 










AP (psl) 

103.3 

12.9 

4. 1 

4. 1 

19.3 

6.4 

4. 1 

12.0 

1 15.8 

4.6 

(kN/m 2 ) 

713. 

09. 

28.2 

28.2 

133. 

44.2 

28.2 

82.8 

799. 

31.7 

E . Port 










K (port)* 

1.6 

1.6 

1.6 

1.95 

. 15 

0.74 

0.50 

1.6 

1.6 

|.6 

AP (ps I ) 

12.4 

3. 1 

1.4 

0.3 

8.2 

0.8 

1.2 

3. 1 

12.4 

0.9 

(kN/m 2 ) 

B5.5 

2 1.4 

9.65 

2.07 

56.6 

5.52 

8.28 

21.4 

85.5 

6.21 

Overa 1 1 









Calculated longest path 
AP (psi) 

199.8 

37.9 

16.5 

12.3 

59.8 

IB. 7 

12.8 

34. | 

205.5 

14. 1 

(kN/m 2 

1378. 

262. 

1 14. 

04.9 

412. 

129. 

88.3 

235. 

14 15. 

97.3 

Kg* Calculated longest path 

479. 

90.6 

39.6 

29.4 

143. 

44.7 

30.6 

78.4 

49 1. 

33.6 

Test effective 

171. 

43.2 

28.5 

19.2 

62.7 

27.3 

20.6 

66.2 

175. 

21.9 


♦See table 10 for areas used in calculating K values 
Air flow = 4 I b/mi n (.0303 kg/sec) 

Density p = .0765 lb/ft 3 (|.225 kg/m 3 ) 









TABLE 10 


LONGEST-PATH PRESSURE DROP ANALYSIS 


Longest path 
components 

Pressure drop 
components 

Flow area based on 
dimensions in table 6 
denoted by item letter 

Inlet Manifold 



A. Port 

90° round bend plus expansion from 
port to fin entrance 

Port area 
n (d) 2 /4 

B. horizontal fin 

Friction for smooth, square ducts 

Horizontal fin flow area 
2 ( i )(g)(b) less fin 
material 

C. Fin miter joint 

Miter bend plus expansion from 
horizontal to vertical fin 

Horizontal fin flow area 
2 ( i )(g)(b) less fin 
material 

D . Ve r l i ca ! fin 

Friction for smooth, square ducts 

Vertical fin flow area 
( i ) (j ) less fin material 

E. Core entrance 

90° miter bend 

Vert 1 cal fin flow area 
( i )(j ) less fin material 

Ovcra 1 1 


Core area 

Outlet Man i fo 1 d 

A . Core ex i i 

90° miter bend 

Vertical fin flow area 
( i )(j ) less fin material 

B. Vertical fin 

Friction for smooth, square ducts 

Vert teal fin flow area 
( i )(j ) less f i n mater ia 1 

C . F i n in t U: r joint 

Miter bend plus contraction from 
vertical to horizontal fin 

Horizontal fin flow area 
2 (i)(g)(b) less fin 
material 

D. Hor i zonta 1 fin 

Friction for smooth, square ducts 

Horizontal fin flow area 
2 ( i )(g)(b) less fin 
ma t e r i a 1 

E. Port 

90° round bend plus contraction 
from fin exit to port 

Port area 
n (d) z /4 

Overall 


Core area 
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Figure I. Thermal Conductivity 



i zoo 1400 1600 1*00 2000 2200 7100 

(667) (77i) (••*) (1000) (Mil) (1222) (1333) 

•turc, °ll (*K) 


Panel Heat Exchanger Materials 
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Figure 2. Transport Properties of Parahydrogen Gas 
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Cooling efficiency, ■ — — = l fl/A)2000 

“actual <" /A) *«“ al 
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Figure 3. Cooling Efficiency 
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Heat flux, q/A 



Figure 4. Heat Exchanger Cross Section Temperature Gradients 






Pressure drop, ab?/i 9 psi/ft (kN/m 3 ) 



Flow rate 


Thermal conductance, T] hA /w l> Btu/sec-°R- f t 2 (kW/°K 







Pressure drop, <rAP/t, psi/ft (kN/m 3 ) 


100 (2040) 


(226) IP 


( 22 . 6 ) 


(2.26) .1 


(.226) .01 


(.0226) .001 
.001 
(.00149 


F i gur 



0 (204) 


(20.4) 


. I (2.04) 
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Thermal conductance, n Q hA T /w£, Btu/sec-°R-f t 2 (kW/°K- 



Pressure drop, aAP/i > psi/ft (kN/m 3 ) 



Figure 8. Effects of Cross Section Shape, Uninterrupted Flow 

Length, and Type of Flow-Length Interruption on Heat 
Exchanger Performance 
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Thermal conductance, T) q h^/wt, B t u/sec- °R- f t (kW/°K 



Thermal conductance, Tl 0 hA T /w|, Btu/sec-°R-f t 2 (kW/ 


(204) 10 


(20.4) 1 


200 (346); 


100 (173) 


IBB 

«■ 

■I 

■I 

II 

II 


3 

H 


ZpjpIO (17.3) 

'30 (51 . 8 ) "~f 


I 


Numbers beside curves are thermal conductivity, 
Btu/hr-°R-ft (W/m-°K ) 

Fin geometry: 

20(7.9)R-. 1 00 (.254)-. 1 25 (.3 1 8 )- .004 ( . 0 1 0 )- ( 10-61) 


.01 

(.0149) 


(.149) 


Flow rate, W/w, Ib/sec-f t (kg/sec-m) 


(1.49; 


Figure 9. Effects of Fin Thermal Conductivity on Thermal Conductance 


8 



No braze fill ets 

.5 in. diam braze 
fillets 


Letters beside curves define fin geometry: 

A. 40TB(l5.8)-.050(. 127)- *4(10. 2)-.005(.OI3)-(7-l,L 0 /4r h =IOO, constant T w ) 

B . 20TB (7.9)-. 1 00 (.2 54)- >8 (20.4 )-.0l0(. 025)- (7-1 , L 0 /4r h = 100, constant T y ) 

C. 40TB ( 1 5.8 )- .050( . 127)- *3 (7 . 6)-.0l0(.025)-(7- 1 , l 0 /4r^= 100, constant T y ) 



Figure 10. Effects of Fin Height, Spacing, and Thickness on Tubular 
Heat Exchanger Performance 
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Thermal conductance, 1 ^ hAy/w£, Btu/sec-°R-f t 2 (kW/ K- 



Pressure drop, <tAP/ 2, ps i /ft (kN/ 


(22 600 ) 1000 


( 2260 ) 100 


( 226 ) 10 


(22.6) I 


( 2 . 26 ) .1 


(. 226 ) .01 



.001 
(. 00149 ) 


.01 .1 
(. 0149 ) (. 149 ) 

Flow rate, W/w, Ib/sec-ft (kg/sec-m) 


( 1 . 49 ) 


Figure II. Rectangular Offset Fin Performance, N = 


15 


£ 


(204 ) 


( 20 . 4 ) 


( 2 . 04 ) 
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Thermal conductance, ^ohA^wi, Btu/sec- 0 R-f t 2 (kW/°K-i 






Pressure drop, ctA P/l, psi/ft (kN/i 


(22 600 ) 1000 


( 2260 ) 100 


( 226 ) 10 


(22.6) I 


( 2 . 26 ) .1 


(. 226 ) .01 



100 ( 2040 ) 


10 ( 204 ) 


I ( 20 . 4 ) 


,1 ( 2 . 04 ) 


.001 

(. 00149 ) (. 0149 ) (. 149 ) ( 1 . 49 ) 

Flow rate, W/w, lb/sec-ft (kg/sec-m) 

Figure 12. Rectangular Offset Fin Performance, N = 20 
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Thermal conductance, T^hAj/wt, Btu/sec-°R- f t (kW/°K-i 





Pressure drop, crAP/l, psi/ft (kN/i 


(22 600) 1000 


(2260) 100 


( 22 . 6 ) 


( 2 . 26 ) . 


(•226) 



E 

i 


10 (204) 


(20.4) 


.001 
(.00149) 


(2.04) 


(1.49) 


Flow rate, W/w, lb/sec-ft (kg/sec-m) 


Figure 13. Rectangular Offset Fin Performance^ N = 25 
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Thermal conductance, T) 0 hA T /wi, Btu/sec-°R-f t 2 (kW/°K 



Pressure drop, aAP/A, psi/ft (kN/m 3 ) 


(22 600) 1000 


(2260) 100 


(226) 10 


(22.6) I 


(2.26). I 


(.226) .01 



10 (204) 


I (20.4) 


. I (2.04) 


.001 

(.00149) 


Figure 14. 


(1.49) 


Flow rate W/w, Ib/sec-ft (kg/sec-m) 
Rectangular Offset Fin Performance, N 


= 30 
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Thermal conductance T^hA^/wZ, Btu/sec-°R-f t 2 (kW/°K-i 





Pressure drop, tAP/£, psi/ft (kN/m 



Figure 15. Rectangular Offset Fin Performance, N = 40 
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Thermal conductance, T1 Q HA /w,^ Btu/sec- 



Pressure drop, oAP, psi (kN/m 2 


Numbers beside curves indicate flow lengths, ft(m) 

Heat exchanger geometry: 20( 7 . 9)R- . 025( . 063) 

- . 05 I ( . 130)-. 004(.0!0)-( 10-59) 



Figure 16. Heat-Flux and Flow-Length Effects on Heat Exchanger 
Performance, h^.. n = 0*025 in. (0.063 cm) 
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Temperature difference, 




Pressure drop, ct AP, psi (kN/m J 


Numbers beside curves indicate flow lengths, ft(m) 



Figure 17. Heat-Flux and Flow-Length Effects on Heat Exchanger 
Performance, h fjn = 0.050 in. (0.127 cm) 
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Temperature difference, 



Pressure drop, crAP, psi (kN/m 2 ) 



Figure 18. Heat-Flux and Flow-Length Effects on Heat Exchanger 

Performance, h.. = 0.100 in. (0.254 cm) 

’ fin 
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Temperature difference, 




Pressure drop, aA P, psi (kN/m 2 ) 



Temperature difference 






Pressure drop, ps i (kN/m 2 ) 




cc 

o 


i 



Heat flux, q/A, Btu/sec-ft 2 (kW/m 2 ) 


Figure 20. Heat-Flux and Flow-Length Effects on Tubular Heat 
Exchanger Performance 
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Temperature difference 



Thermal conductance, T^h/^/wt, B tu/sec- °R- f t (kW/°K 



Figure 21. Effect of Hydrogen Temperature on Thermal Conductance 
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Pressure drop, crAP/A, psi/ft (kN/m 3 ) 

























Legend: 

Valve and hot-temperature sensor 

Low-temperature coolant 

► Intermediate-temperature coolant 

— Maxi mum-temperature coolant 



Structure temperature 



(a) Periodic coolant injection and removal 



(b) Stepwise increment of fin height 

Heat exchanger 



(c) Periodic injection of low-temperature coolant 


Figure 24. Multiple Inlet Concepts 


M- * Valve and hot-temperature sensor 

Low-temperature coolant 
Intermediate-temperature coolant 
Maximum-temperature coolant 



■i | ''-T 'tAt “I 

J, ® 




- h - ,Sr. 

L ♦ 


I) 


Inlet to 
next section 

Outlet 


(a) Two fin rows plus manifold 



Coolant temperature profile 



(b) Three fin rows 


Figure 25. Folded Flow with Injection 


Temperature 


1. Fin geometry: 

16(6. 3)R-. )53(.389)-p1aln-.004(. IOO)-(7-3, L Q /4r h - 100, constant T w ) 

2. W/w = .00934 Ib/sec ft (.0139 kg/sec-m) 

2 

3. o6P < *04 psi (.276 kN/m ) for entire combined length in each series 

4. i nd i cates hydrogen temperature 


— — indicates wall temperatures, and T p 





Figure 26. Temperature Profiles for Fo 1 ded-i n-W i dth Panels of Un 
Unequal Length with Nominal q/A = 10 Btu/sec-ft z 
(114 kW/m 2 ) 
























Design maximum wall temperature; 


Fin geometry 


Btu/sec ft 


conf i gurat ion 


1800 (1000)1 


1700 (944 1 


20(7.9)R-h f . n -l 

L 004 (.0I0)-(| 0-61 ) 
0 

20(7.9)R-h fjn -| 

plain-. 004(.0!0)-(7-3, L Q /4r h ■ constant 

20(7.9)R-h f[n -l 

. 006( . 01 5)- (l 0-61 ) 

0 

20(7.9)R-h f , n -| 

:>lain-.006(.0!5)-(7.3, L Q /4r h « * constant T ) 


500 (5680) 


500 (5680 


250 (2840) 


250 (2840) 


1600 (888) 


1500 (832) 


1. In. ( cm ) 1* n °ted beside points | 

2. Nominal q/A, Btu/sec-ft l (kW/m 2 )i* noted beside 

curves | 

3 - indicates folded and single pass in series 

I— — — indicates two single pass in parallel 


1400 (777) 


Hyd rogen 

parallel 


igure 28 . T r 

























curves (a) and ( c): 

q/A = 250 Btu/sec-ft 2 (2840 kW/m 2 ) 
N =20 fins/In. (7.9 fins/cm) 

1300 (8970)i 1 1 1— | 1 

1200 (8280) 

I 100 (7590) 

1000 (6900) \ 

\ 

900 (6210) ^ 

800 (5520) \ 

700 (4830) r \ l 

600 (4 140) 

| 500 (3450) 

® 400 (2760) 

Q. 

a, 300 (2070) I 1 I 1 I I J 

-Q 

5 (a) Folded with 

- t fln = .004 in. (.010 cm) 


curves (b) and (d); 

q/A = 500 Btu/sec-ft 2 (5680 kW/m 2 ) 
N =20 fins/in. (7.9 fins/cm) 


1300 (978) 
1250 (950) 
1200 (923) 


( b) Folded wi th 

tfj n = .006 in. (.015 cm) 


I I 100 (7590) 

I 

j 1000 (6900) 

I 900 (6210) 
800 (5520) 
700 (4830) 
600 (4140) 
500 (3450) 
400 (2760) 
300 (2070) 



0.020 .060 .100 
(.051 ) (.152 ) (.254 ) 

Fin height, 


( c) Single pass wi th 


004 in. (.010 cm) 


0 .020 .060 .100 
(.051 ) (.152) (.254) 

’fin, in * (cm) 

( d) Si ngl e pass wi th 


1650 (1172) 
1600 (1144) 
1550 (1115) 
1500 (1089) 
1450 (1061) 
1400 (1033) 
1350 (1005) 
1300 (978) 


004 in. (.010 cm ) 


Figure 30. Strength Comparison of High-Flux Folded and Single-Pass 
Heat Exchangers per Figure 27 


Design maximum wall temperature 



(0001 ) 0081 
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Figure 31. Temperature Profiles for Several Flow Routing Concepts 














Temperature difference, 
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( 1000 ) 


Figure 34. Flow Distribution Effects 





(b) Flat tapered 



(c ) Semi cyl i ndrical 



Figure 35. Manifold Concepts 
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Combined inlet and outlet manifold pressure drop, psi (kN/m 2 ) 



(113.5) (1135) (I I 350) (I 13 500) 

(q/A) (i/w), Btu/sec-ft 2 (kW/m 2 ) 

Figure 36. Manifold Pressure Drop 
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Notes : 


I . Full manifold width, w, in. (cm) noted 


beside curves 

2. Port diameter = .8 in. (2.03 cm) 



Figure 39. Calculated Flow Distributions 




Figure 40. 


Calculated Inlet plus Outlet Manifold Pressure Drop 










Core-to-overal 1 total pressure drop ratio 



Core W /W 

max min 


Figure 42 


Effects of Pressure Drop Ratio on Flow Uniformity 


Offset length, 


Numbers beside curves indicate fins per inch and thickness, in. 
(fins/cm, cm) 

Hydraulic radius, r h = .125 (l q ) 


(.406) .16 


(.355) .14 


3 (.3045) .12 


(.254 ) .10 


(.203) .08 


(. 1521 ) .06 


(.1016) .04 



.025 
( . 064 ) 


Fin height, in. (cm) 


Figure 43. Geometry for Rectangular Offset Fins 
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Design maximum wall temperature, 


Numbers beside solid curves indicate flow rates, 
W/w, Ib/sec-ft (kg/sec-m) 


Numbers beside dashed curves indicate hydrogen 
outlet temperatures, °R (°K) 

Fin geometry: 

20 (7 . 9 )- . 075 ( . 191 )-p I ai n- .003 (.008 )- (7-3, 

L 0 /4r = 100, constant T ) 
h W 



T dmw , f , and W/w as Functions of Panel Length for Folded- 
in-Width Panel, at Normal q/A = 10 Btu/sec ft 2 (l|4 kW/m 2 ) 
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Figure 45. 


Intermediate hydrogen temperature, 


N and h r . noted beside curves 
f i n 

Fin geometry: 

(N)R-h.. -plain-. 003(.008)-(7-3, L 0 /4r. = 100, constant T ) 
tin n 



I 18 


Flow rate, W/w, lb/sec-ft (kg/sec-i 



2000 (l | |0 


1800 (1000 


1600 ( 888 ; 


1400 ( 777 ) 


1200 ( 666 ) 


o 1000 (555) 


800 ( 444 ) 


600 ( 333 ) 


400 ( 222 ) 


200 (III) 


Figure 47 . 











Intermediate hydrogen temperature, T , °R (°K) 
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20(7. 9)R-. 1 00 ( . 254 )J 


Folded In width 
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(.0045) 
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(.006) 


C®5) 


,.006 . 
(.0089) 


. 00 7 . 
(.0104) 


,.008 , 
(.01 19) 


Flow rate. W/w, Ib/sec-ft (kg/sec-m) 


Figure 48. T between Fo 1 ded- i n-W i dth and Single-Pass Panels in 

Series at Nominal q/A = 10 Btu/sec-ft 2 (114 kW/m 2 ) and 
A = 2 ft (.61 m) 
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in. (cm) is noted beside curves 


'fin' In. (cm) noted beside curves for nominal q/A of 500 Btu/sec-ft 2 (5680 kW/m 2 ) 


Fin geometry: 


q/A , 

Btu/ sec- ft 
(kW/m* ) 

Flow 

configuration 


250(2640) 

single 

20(7.9 )A-h -L 0 - .004 (.010)- ( 10-6 1 ) 

250(2340) 

folded 

20(7.9)R-h f Jn -plaIn-.004(.0lO)-(7-3 # L 0 /4r h = 100, constant T tf ) 

500 ( 5680) 

single 

20(7.9)R-h f|n -L o -t fln -(lO-6l) 

500(5600) 

folded 

20(7.9)R-h f|n -plaIn-t fin -(7-3, L,/4r h = 100, constant T u ) 



Flow rote, W/w Ib/sec-ft (kg/sec-m) 


Figure 49. T\.. between Fo 1 ded- i n -W i d th and Single-Pass Panels at 

CM 

Nominal q/A - 250 and 500 Btu/sec-ft 2 (2840 and 5680 
kW/m 2 ) and & = 2 ft (.61 m) 









Design maximum wail temperature. 



Figure 50. Fin Geometry Selection for Folded-Flow Only at Nominal 
q/A = 10 Btu/sec-ft 2 (114 kW/m 2 ) and £ = 2 ft (.61 m) 


I 22 
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Figure 51. Operating Conditions and Configurations for Regenerat i ve 1 y 
Cooled Panel Applications 




Letters beside curves designate: 


T co 

W/w 

Allowable hydrogen pressure only 
for fin 20(7. 9)R-.025(. 063)-. 053 
-(. 135) -,003(.008)-( 10-61 ) 

0 R (°K) 

Ib/sec-ft 

(kg/sec-m) 

pslg 

(kH/m* ) 

A. (900 (1058) 

.00314 

(.00460) 

332 

(2290) 

B . 1760 ( 979) 

,00338 

(.00504) 

587 

(4050) 

C. 1600 ( 090) 

.00374 

(.005575) 

852 

(5880) 

0. UOO ( 770) 

,00435 

(.00649) 

960 

(6620) 


Fin geometry (N)r-h - L - .003( . 008)-( 10-6 1 ) 
fin ° 



Fiaure 52. Heat Exchanger Performance for Concepts I and 2 , i - 
9 (.61 m), T r = q/A = 10 Btu/sec-ft 2 (114 kW/m 2 ) 
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2 ft 




Design maximum wall temperature 


Numbers beside curves indicate fins/in. (fins/cm) 
Fin geometry: ( N) R-h^. . -L -• 003( . 008) -( I 0-61 ) 

W/w = .0374 1 b/sec-f t (.0557 kg/sec-m) 



Fin height, in. (cm) 

Figure 53. Heat Exchanger Performance for Concepts 2 and 3, & = 2 ft 
(0.61 m), T r = =°, q /A = 1 00 Btu/sec-ft 2 ( I 1 35 kW/m 2 ) , 

T cq = I 600°R (888°K) 
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Design maximum wall temperature, °F (°K) 


Numbers beside curves indicate fins/in. (fins/cm) 

Fin geometry: (N)R-h -L 0 - .003(.008)-(l0-61 ) 

fin 

W/w = ,1088 lb/sec-ft (.162 kg/sec-m) 



Figure 54. 


Heat Exchanger Performance for Concepts 2 and 3, l - 2 ft 
(0.61 m), T r = q/A = 250 Btu/sec-ft 2 (2840 kW/m z ), 

T = I 400° R (778°K) 

L v/ 
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Fin allowable pressure 





Design maximum wall temperature 



Fin height, in. (cm) 

Figure 55. Heat Exchanger Performance for Concepts 2 and 3, l = 2 ft 
(0.61 m ) } T r = q /A = 250 Btu/sec-ft 2 (2840 kW/m 2 ), 

T CQ = I600°R (888°K) 
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Hydrogen pressure, psi (kN/m 2 ) 
Fin allowable pressure 



Design maximum wall temperature, 



Hydrogen inlet pressure, p s i (kN/m : 





Design maximum wall temperature 


W/w is noted beside curves, Ib/sec-ft (kg/sec-m) 

Tube geometry: 

(N)TB-h f in -pl ai n- . 0 1 0 ( . 025 ) - (7- 1 , L 0 /4r h = 100, constant T ) 



Tube o . d . , in. (cm) 


Figure 57. Tubular Heat Exchanger Performance for Concept 2, 4 = 2 ft 
(.61 m) , T r =«, q/A - 250 Btu/sec-ft 2 (2840 kW/m 2 ) 
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Hydrogen inlet pressure, psi(kN/m 2 ) 
Fin allowable pressure 
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Design maximum wall temperature 



Hydrogen inlet pressure, psi(kN/i 
Fin allowable pressure 








Design maximum wall temperature 



Hydrogen inlet pressure, psi (kN/m 
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Out 1 et temperature, 



Temperature difference 






Outlet temperature 



Temperoture difference,, 








Outlet temperature 


Is noted beside curves 


h 


f i n 


Fin geometry: 


20(7.9)R-h fin -L o -.O03(.0O8)-(lO-6l ) 



300 { 167) 


| 250 (139) 


'200 (III) 


| 150 (83.3) 


100 (55.5) 


■-50 (27.8) 



i 

2 

O 


Flow rate, W/w, Ib/sec-ft (kg/sec-m) 

Figure 63. T co and T ^ - T p for T R = 5000°R (2780°K), q/A = 250 

Btu/sec-ft 2 (2840 kW/m 2 ), Various Flow Rates and Fin 
Geometries with 20 Fins/in. (7.9 Fins/cm) 
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Temperature difference 




Temperature difference^ 






Oesign maximum wall temperature 


N is noted beside curves 

T cq , °R (°K), and W/w, Ib/sec-ft (kg/sec-m), are 
noted above curves 

Fin geomet ry : N (R )-h . - L 0 - .003( ,008 )- ( I 0-6 I ) 
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20 (7.9) 


1760 (976) 


.00296 (.00441) 
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800 (700) | 1 I I I I 

.02 .03 .04 .05 .06 .07 

(.051) (.076) (.102) (.127) (.152) (.178) 

Fin height, h,. , in. (cm) 
fin 

Figure 65. Heat Exchanger Performance for T R = 5000°R (2780°K), 
10 Btu/sec-ft 2 ( I 1 4 kW/m 2 ) 


Temperature difference, 


Fin geometry: 

Concept I, 20(7.9 ) R-. 075 (. 1 90 )- (plain)-. 003 (.008)- (7- 2, L 0 /4r h = 1 00, constant y 

Concept 2, 20(7.9 )R-.025(. 063)-. 060(. 153)-. 003( .008 )-( 10-61 ) 

T =oc. T = 1 600°R (888°K), l = 2 ft (.61 m) 
q > rn 



Figure 70. Concepts I and 2 Off-Design T DMW - T p and aAP/£ for Design 
q/A = 10 Btu/sec-ft 2 ( I l 4 kW/m 2 ) 
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Hydrogen pressure drop, aAP/4, psi/ft (kN/m 3 ) 







Design maximum wall temperature 


F i n geome try: 

Concept I, 20(7.9 )r-. 075(. 191 )- (plain)-. 003 ( . 008 )- (7-2 , L 0 /4r h = 100, constant T ) 
Concept 2, 20(7 . 9 )r- .025( . 063 )- . 060 ( . I 53 )- . 003 ( . 008 )- ( 1 0-61 ). 

T ft oo , T C(J i 600° R (88B°K), J 2 ft (.61 m) 


At design heat flux 


Concept 

P c j » P? 1 a (kN/m* ) 

1 

300.3 (2072 ) 

1 2 

302 *2084 } 



Figure 71. Concepts I and 2 Off-Design T for Design q/A = 10 

DMw ~ 

Btu/sec-ft 2 ( I 1 4 kW/m 2 ) 
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Temperature difference, T ,-T , R (°K) 


2, 40( 1 5 . 7 )R- .038 
Concept 3, 40( t 5. 7 )R- .039 
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• 055 ( . l40)-.003(.008 )-( i 0-6 1 ) 
.056(. !42)-.003(.008)-(|0-6l ) 
= 2 ft (.61 m) 
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Hydrogen pressure drop, oAP/ l 9 psi/ft (kN/i 









Design maximum wall temperature 


Fin geometry: 

Concept 2, 40( 1 5 . 7 )R- .038 (.097 )-.055( . 1 40 )-.003( .008 )- ( 1 0-6 1 ) 
Concept 3, 40( 1 5. 7 )R- .039 ( .099 )-.056(. 1 42 )-.003( . 008 )-( 1 0-6 1 ) 
T r = m > T C0 = |600 ° R (888°K), t = 2 ft (.61 m) 



Figure 75. Concepts 2 and 3 Off-Design T and P„ T for Desiqn 

DMW Cl 

q/A = 500 Btu/sec-ft 2 (5680 kW/m 2 ) 
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Hydrogen inlet pressure, psi(kN/m 2 
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Figure 76. Pressure Drop Sources in Rectangular Manifolds 









Inlet or outlet manifold pressure drop, aA P, psi (kN/m 
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Inlet or outlet manifold pressure drop, oA P, psl (kN/m 1 ) 


(q/A) ( 1 / w), Btu/sec-ft* (kW/m*) Is noted 
beside curves 


Fin geometry: 

10(3. 94) R 250( . 635) -plain-. 004( .010) 
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Figure 78. Density-Adjusted Manifold Pressure Drop vs Manifold Width 
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Port diameter, in. (cm) 


Combined inlet end outlet manifold pressure drop, p$i (kN/m 2 ) 


(q/A) ( £/w), Btu/sec-f t 2 (kW/m 2 ) and in, (cm) 

are noted beside curves 
Fin geometry: 

10(3.94 )R-h fin -piain-. 004 (.010) 

a r = .939 at 500 psia (3450 kN/m 2 ) and |00°R (55.5°K) 
% - .0293 at 250 psia ( 1 725 kN/m 2 ) and |600°R (888°K) 



Figure 79. Manifold Pressure Drop vs Manifold Width 
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Port diameter, in 






Figure 81. Test Setup, Showing Instrumentation and Connected Test Unit 
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Orifice stat i c 
pressure 


Flow measur i ng s 


2 in. 

(5.03 cm) — 

Inlet plenum 
5 static t aps - 


Test unit / 

65 tubes 
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.125 in. (.318 cm) o.d. 
.150 in. (.381 cm) 
center spacing 


Scann i va 1 v 


Bourdon tube 
pressure gage- 


S t at i ons Pi, Pz , P 


Figure 82. Flow Di 








Figure 83. Manifold Test Specimens No. II, 10, 2, and I, Rectangular 
and Tapered Manifold Comparison 
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.2 .4 .6 .8 

Dimensionless core flow width, X/w 


Calibration with Plenum and Core Tubes Only 


Figure 86. 










Static gage pressure, ps i (kN/m 2 ) 



Figure 87. Pressure and Pressure Drop Profiles for Specimen No. 1 
at Inlet and Outlet 


159 


Pressure drop, AP, in. H 2 0 (kN/m 2 ) 





local avg 



igure 90. Flow Profiles for Specimens No. I and II, Rectangular 
and Tapered with One Port Each 



Figure 91. 


Flow Profiles for Specimens No. 2 and 10, Rectangular 
and Tapered with Two Ports Each 



local avg 




Dimensionless core flow width, X/.5w 


Figure 92. Flow Profiles for Specimen No. I, Inlet or Outlet Only 



Core 


IB inlet w 
IA outlet 

i th 

mmm 




m 



mmm 


No. IA inlet 
wi th No. IB 
rout 1 et I 


mmmm 
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Figure 93. Flow Profiles for Specimen No. I with Flow Distribution 
Insert (without Tube Orifices) 
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Figure 94. Longest-Path Manifold Loss Coefficients without Friction 
vs Hor i zon ta 1 - to-Ver t i ca 1 Fin Flow Area Ratio 
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Test effective loss coefficients 



Figure 95. Manifold Pressure Drop Loss Coefficients vs Flow Ratio 
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Figure 97 


Insert Pressure Drop for Uniform Flow in Specimen No 


Core center ] i ne 


Orifice pressure drop, aAP, psi (kN/m 2 ) 





